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Ein grofiziigiges Angebot / A Generous Offer 

Lieber Herr Bohm, 

meine DERIVE-DdXziQn (Version 5.06) zum Unterricht in der Sekundarstufe II (Fach- und Berufsoberschule 
Bayern) habe ich zusammen mit einer kurzen Inhaltsangabe auf meiner neulich erstellten Internetseite in kom- 
primierter Form zum Download bereitgestellt 

http : //lzkopp . gmxhome . de/ index . htm 

Zum Inhalt der Homepage: 

> 35 Dateien in analysis.zip (Analysis fur die 11. und 12. Klasse); 

> 32 Dateien in geometrie.zip (Vektorgeometrie 12. Klasse, Vektoren in Spaltenschreibweise, 3D- 
Grafik); 

> 22 Dateien in stochastik.zip (zB. grafische Simulation von Zufallsexperimenten); 

> 13 Dateien in physik.zip (zB. Auswertung von Messreihen); 

> 38 Dateien in diversesl.zip und diverses2.zip (Was einem so einfallt, nicht nur fur den Unterricht); 

> 3 Dateien in ergaenzung.zip (Radioaktiver Zerfall, harm. Schwingungen, Beugung, Interferenz von 
Wellen). 

Sicher kann mancher DERIVER unter den Fachkolleginnen und -kollegen mit der einen oder anderen Datei 
etwas anfangen. Ich bitte Sie daher, meine Internetadresse im DNL zu veroffentlichen. Vielen Dank! 

Schone Osterfeiertage wunscht Ihnen und Ihrer Familie 
Lorenz Kopp 

DUG-member Lorenz Kopp offers a huge bundle of DERIVE files for download from his website 

http://lzkopp.gmxhonne.de/index.htm 

You can find more than 140 files from Calculus, Geometry, Probability theory, Physics and others. 

I browsed in Lorenz' collection and I was overwhelmed by the plenty of ideas. Many thanks for your 
generous offer. 


GeoGebra - Dynamic Geometry, Algebra and Calculus 

GeoGebra is an open source software that joins 
geometry, algebra and calculus in a new way. It is 
developed for education in secondary schools and 
already received several educational software 
awards. 

On the one hand, GeoGebra is a dynamic geometry system: you can do constructions with the 
mouse and change them dynamically afterwards. On the other hand, you may also work with equa- 
tions and coordinates or calculate derivatives and integrals of functions. These two views are charac- 
teristic for GeoGebra: an expression in the algebra window corresponds to an object in the geometry 
window and vice versa. 

Markus Hohenwarter, University of Salzburg 

www.qeoqebra.at 


GeoGebra 

Markus Hohenwarter 


(GeoGebra is written in many languages. One of its outstanding features is that you can pro- 
duce interactive websites. It is really worth more than only a try. You can download GeoGe- 
bra for free from the above website but you can also work with GeoGebra online. Josef) 
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LETTER Of T Hr rVITGK 


Dear DUG members, 

first of all I'd like to excuse for the delay in deliver- 
ing DNL#57. This is due to two reasons: (1) Easter 
holidays and some courses and seminars which 
needed a lot of preparations and (2) the huge con- 
tents of this DNL. Collecting material for DNL#57 
I intended to add Riideger Baumann's great paper 
on ENGEL-Sequences and Shift registers but when 
inspecting the many mails which I received during 
the last three months I found them so interesting 
that I had to postpone this paper for DNL#58. I 
hope that you will find the questions and answers 
as useful as I do. In my opinion the DNL should be 
a plattform to exchange such findings. 

Let me especially point to the exciting discussion 
on "Stepwise Simplification" (pp 36) and on Albert 
Rich's proposal (p 38). You all are strongly invited 
to participate in improving DERIVE. 

Before adding some notes on this DNL I'd like to 
inform you that you also can download revised 
DNL#5 which is another reason for the delay. It is 
great to reinvestigate the great ideas from 1992 
using means of 2005. Among others there was a 
3D-problem with a poor graphic representation - 
because it was not possible to have a better one in 
1992 - and now we can use implicit and 3D- 
plotting. Some files must be rewritten to be com- 
patible with DERIVE 5 and 6. You can download 
all the stuff. 

What concerns DNL#57 you will find some contri- 
butions removed from the "Preview"-section on 
page 2 - they are published now - but you will also 
find some new contributions announced. I am very 
happy to have a pure Tl-article in this DNL - con- 
tributed by a then 11 year old girl, Mandala von 
Westenholz. At this occasion I'd like to stimulate 
our Tl-users to submit articles. 


Donald McPhillips, busy as ever, provides a paper 
on Actuarial Math ( Versicherungsmathematik ). It 
reminded me that actuarial mathematics was part of 
the curriculum when I was student. One the endex- 
amination tasks in secondary school was from this 
interesting application of probability theory. It 
might be that one or the other would be inspired to 
inform about this applied mathematics chapter. 

Duncan McDougall from Canada joined the User- 
group and sent immediately after transmitting the 
application form an extended paper on Diophantine 
Polynomials. Many thanks. 

finally I'd like to emphasize three highlights of this 
DNL which I strongly recommend to use. On the 
information page you can find two great offers: 

• Lorenz Kopp opens his DERIVE Box and 
invites us all to help ourselves, 

• Markus Hohenwarter offers a great dynamic 
geometry program for free for you and your 
students, 

and on page 49 is the third link 

• Walter Wegscheider worked and is still work- 
ing on Online-HTML-courses for Voyage200 
and DERIVE. 

feedback is greatly appreciated. 

I wish you a wonderful spring time (in the northern 
hemisphere) and a wonderful fall (south of the 
equator) and hope that there will be no delay in 
June for publishing DNL#58 (and revised DNL#6). 

Best regards as ever 



Download all DNL-DERIVE- and Tl-files from 

http : / / www . austromath . ac . at / dug/ 
http : / / www . bk- teachware . com/ main . asp?session=37505 9 
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The DERIVE-NEWSLETTER is the Bulle- 
tin of the DERIVE & CAS- 77 User Group. 
It is published at least four times a year 
with a contents of 44 pages minimum. The 
goals of the DNL are to enable the ex- 
change of experiences made with DERIVE 
and the 77-89 1921 Titanium! Voyage 200 as 
well as to create a group to discuss the 
possibilities of new methodical and didac- 
tical manners in teaching mathematics. 

As many of the DERIVE Users are also 
using the CAS-TIs the DNL tries to com- 
bine the applications of these modern tech- 
nologies. 


Editor: Mag. Josef Bohm 
A-3042 Wtirmla 
DXust 1 
Austria 

Phone/FAX: 43-(0)2275/8207 
e-mail: nojo.boehm@pgv.at 


Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English 
articles will be warmly welcomed nonethe- 
less. Your contributions will be edited but 
not assessed. By submitting articles the 
author gives his consent for reprinting it in 
the DNL. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE & CAS-77 Newsletter will be. 


Next issue: June 2005 

Deadline 15 May 2005 


Preview: Contributions waiting to be published (selection) 

Pringles, B. Grabinger, GER 

Two Stage Least Squares, M. R. Phillips, USA 

Some simulations of Random Experiments, J. Bohm, AUT & L. Kopp, GER 

Wonderful World of Pedal Curves, J. Bohm 

Another Task for End Examination, J. Lechner, AUT 

Tools for 3D-Problems, P. Luke-Rosendahl, GER 

ANOVA with DERIVE & Tl, M. R. Phillips, USA 

Hill-Encription, J. Bohm 

CAD-Design with DERIVE and the Tl, J. Bohm 
Farey Sequences on the Tl, M. Lesmes-Acosta, COL 
Simulating a Graphing Calculator in DERIVE, J. Bohm, AUT 
Henon & Co, J. Bohm 

Challenges from Fermat, Bj. Felsager, DEN 
Are all Bodies falling equally fast, J. Lechner, AUT 
Modelling Traffic Density, Th. Himmelbauer, AUT 

Do you know this? Cabri & CAS on PC and Handheld, W. Wegscheider, AUT 
An Interesting Problem with a Triangle, P. Luke-Rosendahl, GER 
Diophantine Polynomials, Duncan McDougall, CAN 
Rosettes, J. Lechner, AUT 

and Setif, FRA; Vermeylen, BEL; Leinbach, USA; Koller, AUT; Baumann, GER; 
Keunecke, GER and others 

Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wtirmla, DXust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag. Josef Bohm 
Herstellung: Selbstverlag 
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Bill Wilburn WDWNORWALK@aol.com 

ROTATION MATRIX_ANGLES 
HELLO JOSEPH: 

I'M HOPING THAT YOU CAN HELP ME WITH THIS 3x3 ROTATION MATRIX PROBLEM: 

10.3535533905 -0.5732233047 0.7391989197 | 

10.6123724356 0.7391989197 0.2803300858 | 

1-0.7071067811 0.3535533905 0.6123724356 |. 

IT'S A CONCATENATED ROTATION MATRIX, REPRESENTING A ROTATION OF 30 DEGREES 
CCW (1) ABOUT THE X AXIS; 45 DEGREES CCW ABOUT THE Y AXIS, AND 60 DEGREES CCW 
ABOUT THE Z AXIS. 

IF I DIDN'T KNOW THE ROTATION ANGLES, AND I HAD ONLY THE 3x3 ROTATION MATRIX 
GIVEN ABOVE, HOW WOULD I CALCULATE THE ROTATION ANGLES ABOUT THE X,Y AND Z 
AXES? 

NOTE: THE AXIS ROTATION SETS ARE IN ROWS, NOT COLUMNS. 

I HOPE YOU CAN HELP OUT WITH THIS ONE. 

THANKS, 

BILL WILBURN 

(1) CCW = counter clock wise 


Hi Bill, 

I believe that I can help you. 

Please inspect the attached file. It is clear that working with inverse trig functions might cause am- 
biguous solutions, but I think that my process will work. 

Best regards 
Josef 

This was an interesting problem and I am not quite sure if this is the best way to tackle it. But it is a 
nice opportunity to demonstrate the power of a CAS and it might be a challenge for students, too. 

This is to have the concetenated matrix in your form: 

(ROTATE_X (30*1°) ' * R0TATE_V(4 5 . 1 ° ) ’ . ROTATED (60.1°)’)’ 

' ^2 72 3 76 73 ' 

4 $ 4 $ 4 

7e 76 73 3.72 1 

4 8 4 8 4 

72 72 76 

2 4 4 

' 0.3535533905 -0.5732233047 0.7391989197 ' 

0.6123724356 0.7391989197 0.280330065$ 


-0.7071067811 0.3535533905 0.6123724356 
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I coupare with the generalized fom: 

(ROTATE_X(x) ' . ROTATE_r(y) ' . ROTATE J. (z) ' ) ' 

COS(y) *C05(z) SIN(x) *SIN(y) *COS(z) - COS(x) *SIN(z) COS(x) *SIN(y) *C05(z) + SIN(x) .SIN(z) 

COS(y) .SIN(z) COS (x). COS (z) + SIN(x) .SIN(y) .SIN(z) COS(x) .SIN(y) *SIN(z) - SIN(x) .COS(z) 

- SIN(y) SIN(x).COS(y) COS(x) *COS(y) 

I select three eleiaents containing x, y and z and try to solve the respective systen of trig equations: 

[COS (y). COS (z) =0.3535533905, COS (x) .COS (z) + SIN(x) *SIN(y) *SIN(z) = 0.7391939197, COS(x) *COS(y) 
0.6123724356] 

-j2 ^6 73 76 

COS(y) .COS(z) = , COS(x) .COS(z) + SIN(x) .SIN(y) .SIN(z) = + , COS(x) .COS(y) = 

4 3 4 4 

Set cos(x) = x, cos(y) = y and cos(z) = z and go on solving the systen. 

72 2 2 2 76 73 76 ' 

y.z = , x * z + 7(1 - X )*7(1 - y )*7(1 - z ) = — + . x*y = — 

4 3 4 4. 


SOLVE 


42 4s 

y*z = a x*y = , [x, y] 

4 4 


x = 42 * z a y = 


42 


4 * z J 


Substitute for x and y, Simplify and solve for z: 


(73 • z) • z + 7(1 - (73 • z) ) 


4 - 


f 42 2 46 43 

*4(1 - z ) = + — 

( 4*z ) J S 4 


2 2 2 

42*4(1 - z W( 1 - ) *4(8 * z - 1) 2 4s 42 

+ 42* Z = + 

4 ♦ | z | & 4 


SOLVE 


2 2 2 

42*4(1 - z )*4( 1 - 3 *z )*4(s*z - 1) 2 4s 43 

+ 42 * z = + , z , Real 

4 * | z | S 4 


2*4(2742*42 + 10511) 2*4(2742*42 + 10511) 1 1 

z = v z = v z = v z = 

457 457 2 2 

We obtain four solutions for z and proceed with z = 0.5 by substituting in the expres- 
sions for x and y above: 


- 42 *z a y = 


42 


4*z 


x = 73 * a y = 

2 


72 


This is a solution of our problem. 

Accomplishing all the other solutions 
and forming the rotation matrix shows 
that these matrices differ in signs!! 


73 72 

x = a y = 

2 2 


You can follow the process inspecting 
wilburn.dfw (or wi 1 burn5 . dfw) . Josef 



f ^ ) 


72 1 


r 1 Y 

ARCCOS 

<. 2 j 

, ARCCOS 

. 2 > 

, ARCCOS 

v 2 J. 


ACOS 


r 42 


ACOS 


(42 


2 J 


ACOS 


r 1 

^ 2 


= [30, 45, 60] 


1 


1 


1 


= [30, 45, 60] 
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Leopold Watzinger 

I wrote a program to investigate fractals. All single functions are doing fine, but if nesting IF-, ITERTAE 
and VECTOR commands I am failing. You can find my problems in the attached file. 

Expression #14 shall decide the question of convergence for points from [0,0] to [0.2, 0.2] after 20 it- 
erations z n + 1 = z n 2 + / and return the point in case of convergence and [-2,-2] otherwise. 

VECTOR(IF(ITERATE(IF( | [x , y] | < 2, [x - y , 2-x-y + l] , [-2, -2]), [x, y] , [a, 
b], 20) = [-2, -2], [-2, -2], [a, b]), [a, b] , [0, 0], [0.2, 0.2], [0.1, 0.1]) 

DNL: I hope that I can follow your ideas. 

At first I collect all points of the respective region in a list using regi on( . . ) 

#1: region(xw, yw, sw) := APPEND(VECTOR([a, b] , [a, b] , xw, yw, [sw, sw])) 

#2: region([0, 0], [0.2, 0.2], 0.05) 

0 0 ' 

0.05 0 

0.1 0 

Then I include this function into function f ract( ); 


f 2 2 

#4: fractfxw, yw, sw, depth) = VECTOR(IF(ITERATE(IF( | [x , y] | < 2, Lx - y , 2-x-y + 


l], [-2, -2]), [x, y], v, depth) = [-2, -2], [-2, -2], v), v, region(xw, yw, sw)) 
#5: fract([0, 0], [0.2, 0.2], 0.01, 20) 

f 0 0 ] 


0.01 0 

-2 -2 

0.03 0 


This seems to work properly. It could be useful to select all points 
which satisfy the convergence criterium, i.e., which remain with 
abs(z) < 2 after depth iterations. 


#7: fract2(xw, yw, sw, depth) := SELECT(u t [-2, -2], u, VECT0R(IF( ITERATE (IF( | [x , y] | 


[ 2 2 1 
< 2, Lx - y , 2-x-y + lj , [- 


2, -2]), [x, y], v, depth) = [-2, -2], [-2, -2], 


v), v, region(xw, yw, sw))) 


#8: fract2([0, 0], [0.2, 0.2], 0.01, 20) = 


0.01 

0 

0.03 

0 

0.02 

0.01 

0.03 

0.06 

0.13 

0.11 

0.08 

0.13 


0.01 0.16 




p 
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fract2([0, 0], [0.2, 0.2], 0.01, 10) 

DIM(fract2([0, 0], [0.2, 0.2], 0.01, 10)) 

161 

fract2([0, 0], [0.2, 0.2], 0.01, 15) 

DIM(fract2([0, 0], [0.2, 0.2], 0.01, 15)) 

35 

After 10 iterations remain 161 points, after 15 iterations 35 convergent points are remaining and after 
20 there are 8 points left (expr #8). 

We plot the points leaving this filter after 5, 10, 15 and 20 iterations in different colors (red, blue, yel- 
low and black), which leads to a nice fractal picture: 


::::::::::::::::::::: 








Here again we can apply the powerful VECTOR command to produce all the point in one single step. 
The lists can be selected and plotted in different colors, giving nice fractal pictures. 

bigl := VECTOR(fract2([0, 0], [0.2, 0.2], 0.005, k), k, 5, 20, 5) 
big2 := VECTOR(fract2([-0. 2 , -0.2], [0, 0], 0.005, k), k, 5, 20, 5) 


big3 

:= VECTOR(fract2([-0. 2 , 0] , 

[0, 0.2], 

0.005, k), k, 

5, 20, 5) 

big4 

:= VECTOR(fract2([0, -0.2], 

[0.2, 0], 

0.005, k), k, 

5, 20, 5) 


Plot red 

bigl , big2 , big3 , big4 1 

1 1 1 lj 

Plot blue 

bigl , big2 , big3 , big4 1 

2 2 2 2 J 

Plot yellow 

bigl , big2 , big3 , big4 1 

3 3 3 sj 

Plot black 

bigl , big2 , bigS , big4 1 

4 4 4 4j 
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Yves De Racker rc00028433@versateladsl.be 

Mr Boehm, 

As you know the EXPAND function can be used on an expression to expand an expression when it is 
simplified. Recently I tried it with, 

4 3 2 

3-x - 11- x + 9-x + 5-x - 3 

2 

X - 4* x - 2 

(x + ^6 - 2)-(x - -J(> - 2) 

83-x 35 2 

+ + 3-X + x + 19 

2 2 
X - 4*x - 2 X - 4.x - 2 

Instead of two fractions with a linear factor in the denominator! How is that possible? I’m using Derive 
5.06. I send you also a small Derive file where you can see what I have discovered! Thank you very 
much for your help and I wish you a happy New Year for you and your family! 

DNL: 

' 4 3 2 

3-X - 11-X + 9-X + 5 • X - 3 

EXPAND , Radical, x 

2 

v X - 4-x - 2 

67-76 + 166 166 - 67-76 2 

+ + 3-x + x + 19 

4- (x - 76 - 2) 4- (x +7^-2) 


Factorizing the denominator gives : 

But with the EXPAND function you obtain: 


Ernest Carpenter (South Africa) 


Hi Josef, 


I recently looked at some of the latest DNL’s posted on the DUG website and must compliment all of 
you on excellent newsletters full of interesting and useful contributions! 

I’m not a teacher or affiliated with any educational institution (actually to be honest, I’m not even good 
with maths, although I try my best to make sense of some topics©), but I have a great admiration/ love 
for geometry, technology and of course my TI-92+! Thus I enjoy any program, article that involves 
geometry (especially 3D) and the TI-92+. (goes hand-in-hand with my work as CADD Co-ordinator 
and 3D modeler) 


(I also must admit that I’m not a member... yet!) 


A few years back I read a small paper you wrote on representing 3D objects on the TI-92 and you 
have also sent me the program which I used (and is still using) to visualize some basic concepts in 


3D! 


Anyway, the reason for my email is actually to ask you if it would be possible to get my hands on the 
papers mentioned in the “contributions waiting to be published” section of the last few DNL’s: 


- Tools for 3D-Problems, P. Luke-Rosendahl, GER 

- CAD-Design with DERIVE and the Tl, J. Bohm 


They sound intriguing! 

If not possible, would they be in the March DNL? 

Any other references in connection with Tl & 3D geometry that might be useful? 


Thank you for your time, and I hope you have a prosperous 2005! 
Regards 
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Type Checking, Finite Continued Fractions, and 
the Solution of the Linear Diophantine Equation 

ax+by=c 

Stefan Welke, Bonn, Germany , Spwelke@aol.com 

Abstract 

A continued fraction approximation of the quotient f is used to compute a solution 
of the diophantine equation ax + by = c. 


1 Introduction 

Let me first explain the bulky title. The main purpose of this paper is twofold: (1) To introduce the 
function CON_FRAC_RATIONAL(q) which computes the continued fraction expansion of a rational 
number. In this case the expansion is finite. The algorithm is recursive. Thus an irrational argument 
like V2 or n would yield an infinite recursion and a resulting memory overflow. To avoid this di- 
lemma the type checking function rational_Q(a) is introduced. (2) To present a canonical solution 
of the linear diophantine equation ax + by = c as an application of the continued fraction algorithm. 


2 Type checking with DERIVE 

Some algorithms only work with numerical arguments of a special type, e.g.: natural, integer, ra- 
tional, real or complex numbers. The following DERIVE functions return true, if the argument 
evaluates to the type which is specified in the function name. The "Q" at the end stands for query! 


natural _Q(a) := 

If a = FLOOR (a) a a > 0 
true 
false 
false 

natural_Q(i . 5) = false 

integer_Q(a) := 

If a = FLOOR (a) 
true 
false 
false 

[ integer_Q(-10) , integer_Q(^2)] = [true, false] 

[INTEGER?(-IO), INTEGER?Gf2)] = [true, false] 


There is no Expression Type function in 
DERIVE to recognize natural numbers. 

INTEGER?(a) is implemented since 
DERIVE 5. 

DERIVE distinguishes between 
INTEGER? (a) and INTEGER_TYPE? (a) . 
Please inform yourself in the Online- 
Help! 


rational_Q(a) = 

If i nteger_Q( NUMERATOR (a) ) a integer_Q(DENOWINATOR(a)) 
true 
false 
false 


rational_Q(0. 2S999) , 


rational_Q 


& 


2 J 


[true, false] 


- 42 > 

k 2 ^ 


RATIONAL?(0 . 28999 ) , RATIONAL? 


= [true, false] 


D 
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number_Q(a) = 

If DIME NS ION (VARIABLES (a) ) = 0 
true 
false 
false 

[number_Q(2 + k), NUMBER?(2 + k)] = [false, false] 


There is a difference between Stefan's 
number_Q(a) and DERIVE ' s 
NUMBER? (a) treating expressions which 
do not appear as numbers but really 
are numbers. 


[number_Q(2 + k), NUMBER?(2 + k)] = [true, false] 
NUMBER_TYPE?(2 + k) = true 
real_Q(z) :s 

If number_Q(z) a CQNJ(z) = z 
true 
false 
false 


NUMBER_TYPE? (a) is Derivish. 


Note the behaviour of number_Q in the 
examples! 


real_Q(2. 5) , real_Q 


2.5 


= [true, false] 


^ 1 + 2*L 

[real_Q(x + y) , REAL_TTPE?(x + y)] = [false, true] 


real_Q(2.5), real 




= [true, false] 


[real_Q(x + y) , REAL_TYTE?(x + y)] = [false, true] 


symbol ic_Q(a) := 

If DIM (VARIABLES (a)) > 0 


true 

false 

false 


r 2 2i 

L n um be r_Q (2 + m ), symbol ic_Q(2 + m )J = [false, true] 


The functions number_Q(a) and symbol i c_Q are related by the equation 

number_Q(expr)=not(symbol i c_Q(expr)) 


so one of both is redundant. 


3 Continued fraction expansion of rational numbers 

Consider the following chain of equations, 


127 

~46 


■ 2 + 


35 

46 


■ 2 + - 


46 

35 


- 2 H 


1 ~l — — 
1 ^ 35 


:2 + - 


1 + 


3 + - 


5 + i 


which defines the continued fraction expansion of the rational number ^ . The sequence of inte- 
gers [2, 1, 3, 5, 2] contains all information and is commonly used to describe this expansion. We ob- 
tain this sequence by an interative algorithm: 


a a := FLOOR 


127 

~46 


and a n+1 := 


1 


MOD(fl„,l) 


plO 
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The function CONTINUED_FRACTION(a, n) in the NUMBER. MTH file computes n iterations of this 
algorithm, eg.: 


f 127 


#21 : CONT INUED_F R ACT ION 


■, 7 


= [ 2 , 1 , 3 , 5 , 2 , ?, ?] 


46 


This result shows us that we should have terminated the algorithm after four iterations because 
a 4 = 0 mod 1 . The continued fraction algorithm applies to every real number and we have the follow- 
ing 


Theorem 1 The continued fraction algorithm of a real number q terminates iff q is rational 

For a proof see ([1], p. 135-136). We now replace the iteration in the definition of the 
CONTINUED_FRACTION function by recursion. 

cf_vector(v) = 

If V j l = 

vj 2 

c f _vec tor ( [ 1 /MOD (v 1 1 ) , APPEhD(vi2, [ FLOOR (v| 1) ]) ]) 

7 


c f _vec tor 


' 127 
. 46 


[] 


[2, 1, 3, 5, 2] 


As Stefan Welke 's contribution is from pre-DER/\/E-programming times we can now use a 
short program to make the input more comfortable. This is an example how to improve your 
old DER/\/E-functions. 


cf_vector2(x, v) = 

Prog 

v := [x, []] 

#24: If vil = *> 

RETURN v j 2 

c f _vec tor ( [ 1 /MOD (v j 1 ) , APPEM)(vj2, [ FLOOR (vj 1) ]) ]) 


#25: 


cfvec tor 2 


f 127 ^ 


l 46 J 


= [ 2 , 1 , 


3, 5, 2] 


This function terminates iff a n = 0 mod 1 and we need not specify the number of iterations. But we 
must observe that an irrational number would yield infinite recursion and memory overflow. A com- 
bination with a type check prevents us from falling into this trap: 

c o n_f r ac _r a t i o na 1 (a) = 

If rationa1_Q(a) 

#26: cf_vector2(a) 

,J not rational ,J 
,J not rational ,J 


#27: con_f rac_rational (n) = not rational 

r 17 x 


#23: con_frac rational 


L 14 J 


= [1, 4, 1, 2] 
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4 The diophantine equation ax + by = c 

The main purpose of continued fraction expansions is to give good rational approximations of real 
numbers. Let us look at the following example: 

*29: COMT INUED_F R ACT ION ( n , 7) = [3, 7, 15, 1, 292, 1, 1, 1] 

The rational numbers 

1 22 1 333 

[3, 7] = 3 + — = — « 3 .1429 and [3,7,15] = 3 + = *3.1415 

7 7 7 + ^ 106 

serve as very good approximations of n that were already known by Archimedes. Indeed: y- is the 
best rational approximation with denominator m <7 and is the best rational approximation with 
denominator m < 106. For general details see ([1], p. 151). 

Now let us try to find an integer solution of the equation 

127x-46y = 1. 

This equation is equivalent to 

127 y _ 1 
46 x 4 6y 

and we can see that the quotient y must be a very good rational approximation to ^ . In the 
theory of continued fractions the best possible rational approximation is [2, 1,3, 5] = -ff- , i.e. we drop 
the last place of the continued fraction expansion. 

I 53 And surprisingly: 127 • 21 - 46 • 58 = -1. 

2 + = 

. ^ Thus x = -21 and v = -58. This fact is a conse- 

1 + * 

1 quence of the following general 

3 + 

5 


Theorem 2 Let =f a rational number and gcd(a,b) = 1. Let [a 0 ,a l9 ... 9 a n \ 9 n> 1, be the con- 
tinued fraction expansion of £, and let [a 0 ,a l9 . . .,a w ] = be the (n - \)th convergent. 
Then p n \ and q n _\ satisfy 

a-q n - l -b-p n _ l =(-iy. 

For a proof see ([1], p. 13 1 , Theorem 150). Thus continued fractions serve to compute the abso- 
lute values of solutions of diophantine equations of the form ax + by = c with gcd (a,b) = 1. The solv- 
ability of the general case is covered by the next 

Theorem 3 The linear equation ax + by = c with integer coefficients has infinitely many integer 
solutions iff c is a multiple of gcd (af) and no solutions otherwise. Every integer solu- 
tion has the form 
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C Yl 

(x,v)= — ( x 0 ,y 0 + — 7r - rz (-b,a), n eZ, 


gcd (a,b) gcd (a,b) 

where (x 0 ,yo) is an integer solution of gcd " a) x + y = 1 
The proof follows from ([1], p.21, Theorem 24). 


5 Solving ax + by = c with DERIVE 

The implemetation of the solution via continued fractions requires some technical functions which 
shall be briefly explained. It follows from Theorem 2 that the continued fraction solution only works 
if \a\ * \b\. 

We start with two functions borrowed from the file NUMBER . MTH: 

(This is from NUMBER. MTH connected with Derive for Windows 4!!) 

CONVERGENT JWX (v , k_, a) := 

If k_ = 0 
#31: a 

CONVERGENTJWX(v , k_ - 1, vjk_ + 1/a) 

#32: CONVERGENT(v) := CONVERGENTJWX(v , D IRENS ION (v) , «.) 


cf_approx([a, b]) computes the rational approximation of j given as a vector [~a,b]: 

#33: cf_approx(v) := CONVERGENT (DELETE_ELEMENT (v , DIMENSION (v))) 

The next two functions convert fractions '' to vectors [-a,b] and vectors [x,y] to fractions : 

#34: frac_to_vec(a) := [- DENOMINATOR (a) , NUMERATOR (a)] 


#35: 


127 T 


Fracto vec 


L 46 J 


[-46 5 127] 


1 

#36: vec_to_Frac(v) := 

v 

2 

127 

#37: vec_to_Frac([127 , 46]) = 

46 

aux_sol _1 ( [ a , b ] ) computes a solution vector of | ax + by\ = 1 for the case gcd (a,b) = 1 . 

#38: a ux_so1 _1 ( v ) := Fr ac_to_vec (c F_a ppr ox (con_Fr ac_r a t i ona 1 ( vec_to_Fr ac ( v ) ) ) ) 

#39: aux_so1_l([7 , 18]) = [-5, 2] 


co(eq , vars) extracts the coefficient vector [a,b\ from the rationa equation eq := ax + by = c\ 


co(eq, vars) := 


JdCLHSfeq), ^ 
co(27-x + 5-y = 7, [x 5 y]) = [27 3 5] 


vars ), (KLHS(eq), vars )] 

2 J 


aux_sol_2([a, b]) computes a solution of ax + by = gcd (a,b) as a vector. Note that the factor 
v*aux_sol_l(v) corrects the sign of the solution aux_sol_l(v): 
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v ■ aux_so"l_l (v) 

aux_so"l_2(v) := ■ aux_so"l_l(v) 

*43: 2 

GCD(v) 

*44: aux_so1_2 ( [ 127 , 46]) = [-21, 58] 

1 

*45: aux_so1_2([6 , 18]) = 5 0 

. 6 

aux_sol_3(eq) detects the variables in eq and computes a solution of ax + by = c if c is a mul- 
tiple of gcd (af)\ 

#46: aux_so1_3 (eq) := VARIABLES (eq) = aux_so"l_2(co(eq 5 VARIABLES (eq) ))■ RHS(eq) 

#47: aux_so"l_3(10-x + 12 ■ y = 6) = (x = -3 a y = 3) 

#48: aux_so"l_3(10-x - 17 -y = 21) = (x = -105 a y = -63) 

#49: 10- (-105) - 17 -(-63) = 21 

Finally the equation is tested for the condition of the previous theorem and a solution is computed 
if it exists: 

1 i n_d i opha n t_so1 ve (eq ) := 

If HOD (RHS (eq)/GCD (co(eq , VARIABLES (eq) ) ) ) = 0 
#50: a ux_so1 _3 (eq ) 

7 

7 

The function detects the variables but it is crucial to use the correct order ax + by = c, i.e. vari- 
ables on the left and numbers on the right hand side, and one must observe the conditiot \a\ ^ \b\: 

#51: 1"in_d"iophant_so1ve(- 127 -s + 46«t = 5) = (s = 105 a t = 290) 

#52: 1"in_d"iophant_so1ve(- 12 -x + 18 -y = 15) = ? 

#53: 1nn_dnophant_so"lve (1023797 -u + 4601337 -v = -15078) 

#54: u = -29117563062 a v = 6478654728 

#55: 1"in_d"iophant_so1ve(- 127 -s = - 46«t + 5) = ? 

#56: 1 - in_d - iophant_so1ve(x + y = 3) = (x = -m a y = m) 

Though -127s = -46 1 + 5 is equivalent to -127s + 46^ = 5 this function in its present form cannot 
solve the first equation. 

Let me finally explain the reason that makes the continued fraction solution so particular. 
ax + by = c is the equation of a straight line in [R 2 , which is perpendicular to the vector ( a , b). All in- 
teger solutions are points on this line with integer coordinate values. The following theorem explains 
the geometric significnace of our solution: 

Theorem 4 Let the continued fraction solution (x, y) of ax + by = c exist. Then (x, y) has the least 
distance to the origin among all integer solutions. 

I have no reference for a proof of this theorem but I cannot imagine that no proof in the literature 
on number theory exists. 

We already mentioned the impossibility to construct a continued fraction solution in the case 
\a\ = \b\. See expression #56 above. 

The arithmetic reason for this failure is the following fact: The continued fraction expansion of 
Ifl [1], i-e. we cannot approximate by a shorter continued fraction. In geometric terms this differ- 


ence is: 
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• If \a\ ^ \b\ and an integer solution exists then exactly one integer solution minimizes the 
distance to the origin. 

• If \a\ = \b\ and an integer solution exists then exactly two integer solutions minimize the 
distance to the origin. 

A piece of additional code that chooses one of the two minimizing solutions closes this gap. We 
resemble in this way the geometric flavour of the continued fraction solution. 


c - v -FLOOR 

1 2 ■ v 


aux_so"l_4(v , c) := lx = FLOOR 


#58: aux_so"l_4([l , 1], 3) = [x = 1, y = 2] 

#59: aux_so"l_4([3 , -3], 9) = [x = 1, y = -2] 


The next function tests for the condition \a\ = \b\. 

eq_test(v) := 

IF ABS(vU) = ABS (vj2) 

#60: true 

False 
False 

We now enhance the solving function: 

lin_diophant_solve(eq) := 

If NOD (RHS (eq) /GCD (co(eq , VARIABLES (eq))) s 1) = 0 
If eq_test(co(eq s VARIABLES (eq))) = true 
#61: aux_sol_4(co(eq , VARIABLES (eq) ) s RHS(eq)) 

aux_sol_3(eq) 

7 


#62: lin_diophant_solveO ■x + 3-y = 18) =[x=3,y=3] 

#63: lin_diophant_solveC3 -x + 3-y = 178) = ? 

#64: lin_diophant_solveC3 -x - 3 -y = 9) = [x= 1, y = -2] 

#65: lin_diophant_solve(3-x + 5-y = 178) = (x = 356 a y = -178) 

#66: lin_diophant_solve(6-x + 18-y = 18) = (x = 3 a y = 0) 

#67: lin_diophant_solve(6-x + 6-y = 18) =[x=l s y=2] 

Now the solving function covers all cases. 

Remarks: (a) the method of solving linear equations with continued fractions is briefly described 
in ([2], p.195). The continued fractions solution is restricted by its nature to equations with two vari- 
ables. (b) Historically the theory of continued fractions was further developed by Lagrange to solve 
PelTs equation x - dy 1 = 1 with d e N and yfd € N. This method and its treatment with DERIVE wil 
be presented in a forthcoming paper. 
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This worksheet provides a set of basic functions, such as D, N, S, C, M, and R, for comput- 
ing the actuarial present values of life insurance and life annuities, as well as premiums and 
contributions. It assumes a basic working knowledge of actuarial mathematics and Derive. 

Each actuarial function allows for calculations based on three different life tables: a primary 
life (X), a secondary life (Y), and a joint life (J). Mortality for the primary and secondary lives 
is specified by selecting appropriate mortality tables. A good source of hundreds of mortality 
tables is The Society of Actuaries Table Manager 3.0 sofware that can be found at 
www.soa.org. Simply enter Table Manager in the search box on the home page and follow 
the instructions for downloading and using the software. The individual mortality tables can 
be exported from the Table Manager to a text file. If you cut out all the extraneous text and 
save only the age and mortality numbers in a text file, the file extenstion can be renamed to 
.dat and loaded into Derive as a data file. The resulting 2-column matrix should be given an 
appropriate name and then saved as a .mth file for later use with the functions in this work- 
sheet. 

I used two sources in developing the actuarial functions. The first is Actuarial Mathematics, 
by Bowers, et. al., The Society of Actuaries, 1986. The second is QLIB/48: Users Guide, by 
Steve Lindblad, 1997; this is an actuarial library for the HP-48GX calculator. Any errors in 
the concepts or code developed here are, of course, my own! If, by chance, you find any 
errors, please let me know at one of the email addresses given above. 

On the other hand, if you find this worksheet useful, I would also appreciate hearing from 
you. Any comments, questions, or suggestions for improvement are, of course, most wel- 
come. 


The functions and their arguments included in this worksheet are (where x means a person 
age x): 

1. CommutationFunctions(tabs , i , axmy : =0 , too , r:=10 A 9) , compute commutation tables 

2. lx(age,tab:=x) , number living at age x 

3. Qx(age,tab:=x) , probability of x dying in one year 

4. Px(age,tab:=x) , probability of x living one year 

5. nPx(age,n,tab:=x) , probability of x living n years 

6. nQx(age,n,tab:=x) , probability of x dying within n years 

7. tnQx(age,t,n,tab:=x) , probability of x living t years and then dying within 
n years 

D(age,tab:=x) , D commutation function 


8 . 


pi 6 


MacDonald R. Phillips: Actuarial Math with DERIVE 6 


D 


-N-L#57 


9. DM(age,m,tab:=x) , DM commutation function 

10. N(age,tab:=x) , N commutation function 

11. NM(age,m,tab:=x) , NM commutation function 

12. S(age,tab:=x) , S commutation function 

13. SM(age,m,tab:=x) , SM commutation function 

14. C(age,tab:=x) , C commutation function 

15. CM(age,m,tab:=x) , CM commutation function 

16. M(age,tab:=x) , M commutation function 

17. MM(age,m,tab:=x) , MM commutation functions 

18. R(age,tab:=x) , R commutation function 

19. RM(age,m,tab:=x) , RM commutation function 

20. CaM(n,m), n-year certain annuity payable mthly 

21. CLM(age,n,m,tab:=x) , n-year certain & life annuity payable mthly 

22. CL12(age,n,tab:=x) , n-year certain & life annuity payable monthly 

22. 3SM(age,p,m) , joint & survivor life annuity payable mthly 

23. DS12(age,p), join & survivor life annuity payable monthly 


1. Setting the Assumptions 

The commutation functions are determined by the mortality table(s) you select, any age set- 
backs or set forwards, the interest rate, age difference between the primary and secondary 
lives, terminal age of the table(s), and radix. All this information goes into the Commutation- 
Functions(tabs,i , axmy : =0 , too , r:=10 A 9) program. 

The arguments are: 

tabs -- a vector or 2-row matrix specifying a primary or a primary/seconday mortality table 
and any age setbacks or set forwards. Setbacks or forwards allow you to shift the mortality 
rates up to later ages (setback) or down to earlier ages (set forward). For example, a 2-year 
setback shifts the mortality rates up two years so that a person 50 years old, for example, will 
now have the same mortality rate as a person 48 years under the original, unshifted table. 
Setbacks therefore decrease the level of mortality for any given age (if the mortality rates are 
increasing, as is usually the case). Set forwards, on the other hand, increase the level of 
mortality for any given age. Setbacks are entered as positive numbers and set forwards as 
negative numbers. If there is no setback/forward simply enter a 0. 

i -- the assumed interest rate expressed as a percentage, i.e., an interest rate of .06 is en- 
tered as 6. 

axmy -- the age difference between x and y. For example, if x is 62 years old and y is 60, 
axmy = 2. On the other hand, if x is 60 and y is 62, axmy = -2. The default for axmy is 0 (zero). 

too -- the terminal age truncation options allows you to set an age beyond which the table are 
not permitted to extend. If the terminal age of the mortality table is greater than the too value, 
mortality rates at the end of the table are discarded to force the desired terminal age. This 
can be useful when trying to match commutation functions produced with software based on 
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truncated tables. For example, if too is set to 100, then mortality rates above 99 are dis- 
carded and the mortality rate for age 99 will be set to 1 , so that the terminal age of the table 
is 100. If no value is set for too or it is set to an age greater than the truncation age of the 
table, no truncation will be done. 


r -- the radix of the table. The radix is the initial number of lives assumed to exist at the ini- 
tial age of the mortality table. The default value is 10 A 9. 

The mortality tables used in the examples are the Male UP-94 and Female UP-94 tables. 
"UP" refers to uninsured pensioner. They are tables #832 and #833 in the Society of Actuar- 
ies Table Manager 3.0 software. The tables are named mupl994 and fupl994 in this work- 
sheet. They may be viewed by going to Author/Variable Value... and selecting the name of 
the table. 


These two tables are used with no age setbacks or forwards, an interest rate of 6 percent, an 
age difference of 2, no terminal age truncation, and the default radix. 


# 1 : 


Commu tat ionFunct ions 

mupl994 0 

N 

, 6, 2 

\ 

. fupl994 0 . 

y 


Notice that a 2-row matrix was used to specify the mortality tables and age set- 
backs/forwards. In this instance, there are no setbacks/forwards, they are set to zero (0). In 
specifying two mortality tables, a primary and secondary table of commutation functions are 
computed, as well as a joint-life table calculated from the primary and secondard tables and 
based on the age difference you specify. If you want to use only one table based on the pri- 
mary life, simply enter a 2-element vector speficifying the mortality table and age set- 
back/forward. Any prior secondary and joint tables will be deleted. 


Now, select the function and click on the = (Simplify) button in the toolbar. 

#2: Tables Done! 

Since the tables are very large, they are not displayed. They have, however, the global na- 
mes of x, y, and j. The first four rows of the x (primary) table are displayed below. The first 
column is the age, the second is the number living at the beginning of that age, the third is 
the D commutation function, the fourth is the N commutation functions, the fifth is the C 
commutation function, and the last is the M commutation function. In general, the D and N 
commutation functions are used for life annuities and the C and M commutation functions for 
life insurance. All calculations done with the commutation functions use a ratio of at least 
two of the functions. 


Two other global variables created are i rate -- the interest rate, and qxmy -- the age differ- 
ence between X and Y. 
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X 

#3: [1 4] 


Age 

1 

D 

N 

C 

M 

1 

1000000000 

943396226.4 

16344980654.4 

566927.7322 

18208642.2 

2 

999363000 

889429512.2 

15401584427.9 

360806.3115 

17641714.47 

. 3 

998933273.9 

838723639.2 

14512154915.7 

282475.7917 

17280908.15 . 

irate 







# 6 : 6 
#7 : qxmy 

# 8 : 2 


2. Using the Commutation Functions 

Example 1 : Compute the actuarial present value (APV) of a life annuity of 1 payable at the 
beginning of each year to X age 50 based on the assumptions in the commutation functions 
computed above. (By stipulating "X age 50" I mean a person age 50 using the primary table 
of commutation functions.) 

The APV of a life annuity of 1 payable at the beginning of the year is given by Nx/Dx, i.e., the 
N commutation function for age x divided by the D commutation function for age x. 

N(50) 

#9: 

D(50) 

#10: 13.99577103 

The APV is 13.99577... Notice that in the N and D commutation functions I did not specify a 
primary or secondary life; this is because the commutation functions default to the primary 
life (X). If I wanted to compute the same APV for the secondary life (Y), I would enter: 

N(50 , y) 

# 11 : 

D(50, y) 

#12: 14.78052895 

Or, for an annuity payable on the joint life status of X and Y, i.e., as long as both X and Y are 
alive. 

N(50 , j) 

#13: 

D(50, j) 

#14: 13.30680737 

Example 2: Compute the same APV but with the annual annuity of 1 payable monthly. That 
is, the annuitant receives 1/12 of the annuity every month. For this caculation use the 
NM(age,m,tab:=x) function instead of the N function. In this case, m would equal 12 since 
the annuity is paid monthly, m can take on any value from 1 to 8; 8 is used for annuities 
payable continuously. For X age 50 we have: 
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Example 2: Compute the same APV but with the annual annuity of 1 payable monthly. That 
is, the annuitant receives 1/12 of the annuity every month. For this caculation use the 
NM(age,m,tab:=x) function instead of the N function. In this case, m would equal 12 since 
the annuity is paid monthly, m can take on any value from 1 to 8; 8 is used for annuities 
payable continuously. For X age 50 we have: 

NM(50, 12) 

#15: 

D(50) 


For Y age 50 we have: 


#16: 13 . 53158441 


NM(50, 12, y) 

#17: 

D(50, y) 

#18: 14.31656285 


Example 3: Compute the same APVs but with the annuities payable at the end instead of 
the beginning of the period. For the annuity payable annually simply increase the age of the 
N commutation function by 1 . For X age 50: 

N(51) 

#19: 

D(50) 

#20: 12.99577103 

To calculate the APV of annuities payable mthly at the end of the period, enter m as a nega- 
tive value. 

NM(50, -12) 

#21 

D(50) 

#22: 13.44825108 

(The fact that a negative payment frequency gives you annuities payable at the end of each 
period is a consequence of the mathematical definition of the actuarial functions.) The result 
given by N(5l)/D(50) can be gotten, of course, from setting m to -1 in the NM function. 

NM(50, -l) 

#23: 

D(50) 

#24: 12.99577103 


Example 4: Calculating a 75% Joint-and-Survivor Annuity 

The j table, along with the x and y tables can be used to calculate joint-and-survivor (or last- 
to-die) annuities. Assume X, age 65, receives an annuity of $30,000 a year as long as he 
lives and his wife Y receives 75% of $30,000 (or $22,500) for as long as she lives after X 
dies. Assume also that the annuity is payable monthly at the beginning of the month and 
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that Y is 2 years younger than X. Use the commutation functions already computed with the 
mortality tables, setbacks/forwards, and interest rate as stipulated above. Compute the APV 
of this annuity. 


This APV can be computed in 3 steps. First, X will receive an annuity as long as he lives of 
$30,000. The APV is: 

30000NM(65, 12) 

#25: 

D(65) 


#26: 303285.7073 

Second, Y will receive an annuity of $22,500 as long as she lives. Since she is 2 years y- 
ounger than X, enter her age as 63. The APV is: 

22500NM(63, 12, y) 

#27: 

D(63, y) 

#28: 265181.3315 

But Y will not recieve her annuity as long as both X and Y are alive. From the addition of 
their annuities' APV, we need to subtract the APV of $22,500 for as long as both are alive. 
The joint mortalilty is indexed to X's life with Y's age equal to X's minus their age difference. 
So, when using the joint commutation functions enter only X's age. The APV is: 

22500NM(65, 12, j) 

#29: 

D(65, j) 

#30: 202485.1758 


The APV of this 75% joint-and-survivor annuity is therefore $365,982. 

#31: 303285.7073 + 265181.3315 - 202485.1758 = 365981.863 

If X is alive, he receives his annuity of $30,000; Y's annuity and the joint annuity of $22,500 
cancel each other out. If X dies and Y remains alive, X's annuity of $30,000 is no longer paid 
and the joint annuity of $22,500 is no longer viable either; this leaves Y with an annuity of 
$22,500 for as long as she lives. The same APV can be computed with the DS12(age , p) 
function, where p is the percentage of X's annuity received by Y. 

#32: 30000 -3S12(65, 75) 

#33: 365981.863 


Example 5: Suppose X, age 25, wants to provide for a retirement annuity of $2,500 per 
month when he retires at age 65. What is the APV of this deferred annuity? 

2500-12-NM(65, 12) 

#34: 

D(25) 


#35: 25784.22814 
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If X wanted to make annual payments to fund his deferred annuity, what would they be? 

2500-12-NM(65, 12) 

#36: 

N(25) - N(65) 

#37: 1641.925657 

X's contribution would be $1 ,641 ,93 payable at the beginning of each year until he was age 
64. OrX could make an annual contribution of $1,687.48 with 1/1 2th ($140.62) payable at 
the beginning of each month until he was age 64 and 1 1 months. 

2500-12-NM(65, 12) 

#38: 

NM(25 , 12) - NM(65, 12) 

#39: 1687.477924 


Example 6: Y wants to receive a 10-year certain and life annuity of $3,000 per month pay- 
able at the beginning of the month when she retires at age 65. What is the APV of the annu- 
ity at age 65? 


This is a combination of a 10-year certain annuity and a 10-year deferred life annuity. The 
function CaM(n,m) computes the present value of an annual annuity payable mthly for a cer- 
tain number of years; if m is positive, it is an annuity due. If m is negative, it is an annuity 
immediate. And, if n is negative, CaM(n,m) computes the future value of the anuity certain. 
The APV of the certain and deferred life annuity, payable monthly, is $421,531. 


#40: 


#41: 


3000 -12- 


CaM(10, 12) + 


v 

421530.7227 


NM(75, 12, y) > 
D(65, y) „ 


If instead, the annuity was to be paid at the end of the month, simply enter -12 for m. The 
APV is $418,763. 


#42: 


3000 -12- 


CaM(10, -12) + 


V 


NM(75, -12, y) " 
D(65, y) 


#43: 418762.6292 


Two other functions included in the worksheet will compute a n-year certain and deferred life 
annuity. They are CLM(age,n,m,tab:=x) and CL12(age,n,tab:=x). 


#44: 300012CLMC65, 10, 12, y) 

#45: 421530.7227 

#46: 3000- 12 -0.12(65, 10, y) 


#47 


421530.7227 
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Example 7: A firm's pension plan provides at retirement at age 65 $50 per month per year of 
service for its employees. It uses the Entry Age Normal (EAN) method of funding it's pension 
plan, an interest rate of 6% and the UP-1994 mortality tables. X, age 30, is just hired. What 
is the firm's normal cost, i.e., annual cost, for funding X's expected pension in the year he is 
hired and for all subsequent years assuming no changes in the retirement plan? 


In this case, the normal cost is based on X's expected penion at age 65. This is $50 * 12 
months * 35 years of service = $21 ,000. The normal cost is: 


#48: 


#49: 


nc • (N(35) - N(65)) 


21000-NM(65, 12) 


SOLVE 


D(65) 
nc • (N(35) - N(65)) 
0(65) 


D(65) 

21000NM(65, 12) 
0(65) 


A 

nc , Real 

y 


#50: nc = 2276.983661 


In this example, the APV of the pension at age 65 (right hand side of #48) is equated to the 
actuarial future value at age 65 (left hand side of #48) of the annuity used to fund the pen- 
sion. 

Example 8: What is the APV or net single premium of a whole life insurance policy of 
$100,000 payable at the end of the year of death for X age 25? The APV of a life insurance 
policy is given by Mx/Dx. 

100000M(25) 

#51: 

D(25) 

#52: 6022.617956 


If X wants to pay annual premiums on this policy at the beginning of each year for as long as 
he lives, what would it be? 


Monthly premiums would be: 


100000M(25) 

#53: 

N(25) 

#54: 362.7499466 

100000M(25) 
#55: NM(25 , 12) 

12 


#56: 


31.09722343 
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Example 9: Y wants to insure her life for $1 ,000,000 during her working years from age 25 to 
65, payable at the moment of death. But she only wants to pay annual premiums for 20 
years. What is the amount of the premium? 


For this example use the MM(age , m , tab : =x) function with m set to 8. For life insurance, m 
determines when the benefit is paid. A positive number means it is paid at the end of the 
period of death. A negative number means it is paid at the beginning of the period of death. 
For example, if m = 12, the death benefit is paid at the end of the month of death. If m = 1, it 
is paid at the end of the year of death. If m = 8, it is, of course, assumed to be paid at the 
moment of death. If you noticed above that the sign of m means the reverse for life annui- 
ties, you are correct. Again, this is a consequense of the mathematical definition of the actu- 
arial functions. 


1000000 (MM(25, «, y) - MM(65, «, y)) 

#57: 

N(25 , y) - N(45 , y) 

#58: 1230.800747 

Example 10: The commutation functions S(age,tab:=x) and R(age,tab:=x) are used for 
calculations involving increasing annuities and increasing life insurance, respectively. 

SM(age , m , tab : =x) and RM(age , m , tab : =x) refer, of course to increasing life annuities and life 
insurance payable mthly. For instance, what is the APV of a life annuity of $1 ,000 that in- 
creases by $1 ,000 every year to X who is age 65? 

1000 S(65) 

#59: 

D(65) 

#60: 90139.53494 

Or, what is the APV of a life insurance that pays $10,000 the first year if death occurs and 
increases by $10,000 each year after to Y who is age 35? 

10000 R(35, y) 

#61: 

D(35, y) 

#62: 28970.10445 

3. Non-integral Ages and the Uniform Distribution of Deaths 

All the commuation functions accept non-integral ages, e.g, 50.5 or 32.3. Their values are 
calculated on the assumption that deaths during a year follow a uniform distribution. How- 
ever, if non-integral ages are used with functions such as NM or MM where abs(m) > 1, there 
are small relative errors in the calculations. Towards the end of the mortality table the errors 
can become large. If non-integral ages are used only when abs(m) = 1 ; or if only integral 
ages are used when abs(m) > 1 , there is no error, given the uniform distribution of deaths 
assumption. 
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4. Q's and P's 

The Q's and P's associated with actuarial math are the probabilities of dying and living. (In 
actuarial textbooks they are always lower case. Since DERIVE does not do sub- or super- 
scripts very well, I have capitalized Q and P.) For instance, based on the mortality table, Qx 
is the probability that a person age x will live one year. Px is the probability that a person 
age x will live one year. nQx is the probability that a person age x will die within n years; nPx 
is the probability that a person age x with live n years. And finally, tnQx is the probability that 
a person age x will live t years and then die within n years. (The function lx(age,tab:=x) gives 
the number living at any given age. It is used to calculate the Q's and P's. 

The values of the q's and p's are available for X, Y, and J, that is, the primary life, secondary 
life, and joint life status of X:Y. 


Example 1 1 : What is the probability that X age 50 will die in one year? Ten years? Will live 
10 years and then die in the next 5 years? 


#63: Qx(50) 

#64: 0.002773 

#65: nQx(50, 10) 

#66: 0.04691685322 

#67: tnQx(50 , 10, 5) 

#68: 0.0516509792 


Example 12: What is the probability that the joint & survivor status of X:Y will last 10 years if 
X is 65 and Y 2 years younger? Since qxmy = 2, we have: 

#69: nPx(65, 10) + nPx(65 - qxmy, 10, y) - nPx(65, 10, j) 

#70: 0.97409016 

The probability that either X or Y or both will be alive in 10 years is .974. 

This is simply the probability that X and Y survive 10 years minus the joint probability that 
they both live 10 years. Subtracting the joint probability basically eliminates double counting. 
See any textbook on elementary probability theory. 

5. Conclusion 

This is just a taste of what can be done with actuarial mathematics. The commutation func- 
tions can be used for calculating life insurance premiums, reserves and cash values, or the 
normal cost and actuarial value of pensions, or contributions to tax deferred annuities, etc., 
etc. 
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AVOIDING CONVOLUTION AND TRANSFORM METHODS 

By: Milton Lesmes Acosta. Bogota, Colombia. 

1. Intersection of areas and Convolution 

When you need to express the amount of overlap of one object over another you will find some com- 
mon aspects with the mathematical concept of convolution. 

The following picture illustrates the area of intersection of a static triangle and a square which moves 
in the direction of V. 



Area 

of overlap 



Compare this with the measure of the overlap of two real integrable functions /, g defined over the 

t 

interval [0,/] and the expression J/(x)g(/ 1 -T)<ir 

o 

2. Linear Time Invariant Systems. 

The notion of convolution appears naturally in Linear Time Invariant Systems. 

In the distribution theory the representation 

+oo 

/(»)= J/(t)5(»-T)rfc 


is a tool to explain why the convolution integral appears in this context. The signal / is processed by 
a LTI System H affecting the Dirac delta distribution because of the linearity 
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+00 +00 

= H( { /( X )5 (t -T ]dx) = J /(x )H{ 8 (f -T )) A 


—00 


—00 


Let the response of the System to the Dirac delta distribution be H {p ( t )) = h{t ) then if the LTI 
System is invariant in time (that is the system always responses in the same way no matter the instant 
of time) it follows that H( s(»-t)) = h{t -T ) , then we get 


H(f(t))= \f{x)h{t-x)dx 


which is known as the convolution between / and h 


+oo 


f*h= | /(x )h(t -x )dx 

—00 


3. Convolution and DERIVE 


For a Derive user in many cases it is easy to calculate this integral, for example if 

fit ) := t.chif),t,l) 
git) := 2.c/j/(0,t,2) 

then 

INT(f(x )-g(t -x ), x , -oo, +oo ) 
is the convolution represented in the following picture 



4. The discrete convolution. 

The convolution of two discrete signals is given by 

+oo 

fim)*gim)= Yj f(k)gim-k ) 

k =- oo 
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For two discrete signals /, h defined over the set {l, 2,3,..., «}we have [/(l),/(2),/(3), ...,/(«)] 
and [h(l),h(2),h(3),...,h(n)] the derive user can get the convolution f *h as the coefficients of the 
product 

(/( l)x + /( 2)x 2 + /(3)x 3 +...+ /(/7)x w ^/z(1)x + /2(2)x 2 +/?(3)x 3 + 


Maybe it is important to remember the DNL solutions to Rainer Wonisch's question in Newsletter #49 
p 4. and the findings of some DERIVER's. Of course, you will recognize the z-transform in this proc- 
ess and the interested reader must see the dfw files of Terence Etchells. 

One way to get the coefficient of x m of the polynomial p(x) is: 


define p(x), calculate 


and evaluate at x:=0. 


DIF(p(x),x,m) 

m ! 


5. The sum of two independent random variables and convolution. 

In the random experiment consisting of tossing a fair die, for the random variable X “the number of 
the point in the die”, the probability mass function is shown in the following figure 


The random experiment of tossing two fair dice has a probability mass function which can be calcu- 
lated as a convolution to get 
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This was done with the coefficients of the product 


1 1 2 1 3 1 4 1 5 1 6 

— X H X H X H X H X H X 

\6 6 6 6 6 6 j 


1 1 2 1 3 1 4 1 5 1 6 

— JCH JC H JC H JC H JC H JC 


v6 6 6 6 6 6 

In the continuous case with the uniform distribution that has the probability mass function 

1, for 0 < x < 1 

/>) = 

10, otherwise 

the probability mass function of the sum of the random variables each distributed uniform is the con- 

X 

volution f = f x * f x -> the distribution function F(x) = J / ( t)dt is shown in the following picture 


As it is seen it is needed to calculate a convolution and a integral to get the distribution function. Of 
course that for the integral we could use a Monte Carlo method. 


6. The Derive, avoiding Convolution and Monte Carlo Method. 


First it is described the method to be used: Take samples of the random variable. Sort this samples. 
Distribute this samples over the interval (0,1) 

This is done with the derive instructions 


W := sort(vector ([s(£)] , k, 1, 1 0000) ) 


vector 


det([^]), : 


,z, 1,1 0000 


10000 _ 

Where s(k) is the k-value of sample of the random variable, it is needed to implement a method to 
generate s(k) depending of the random variable (or the experiment). 


Example: Choose two independent numbers of the interval [0, l] calculate the sum (that is the continu- 
ous random variable of example of 5.), The 10000-sample of this random variable is in 
W := s ort (vector ([r andom(l) + random(l)],k, l, 10000)) 
rr 1 \ 


det([^]), 


l 

10000 


,z, 1,1 0000 

) 


is in the following graph 


and the distribution vector 
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m ; r; r 




this is a nice result, and it is seen how the convolution and the counting of the generated points (to 
apply the Monte Carlo method) are avoided. 

The following figure shows a distribution of 10000 points result of the instruction 


vector 


( r k 

\ 

ranaomCl ) + ranaomCl), 

, it, 1,1 0000 

ioooo 

) 



7. Trying with discrete cases. 


Finally it is shown the result of the formula vector 
the random variable is discrete. 


det([^]), 


10000 


,1,1,10000 


in case where 


W := sort (vector {[random(6) + 1 + random{6 ) + 1], £,1,10000)) 
correspond to a 10000-sample of the random experiment of tossing two fair dice, the following figure 


shows the result of applying vector 


det([^]), 


10000 


,2,1,10000 


to this sample. 
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Martin Lindsay MLindsay@distance.vic.edu.au 

Dear Joseph Bohm 

I have attached my DUG membership form. 

I was wondering whether you could help me. My university has just purchased Derive 6. 

My problem is this. I’m doing a first year Calculus unit. The usual stuff that you find in any US 
$130 textbook(l): Derivatives, Applications of, integrals, techniques of, Inverse Trig, Hyper- 
bolics, etc. I thought, given my background and research (I have a PhD in mathematics edu- 
cation) I could put together some Derive labs on this stuff, but it is proving a little harder than 
I thought. 

For a start, there seem to be no Derive/CAS books as such at this level, most of it is pre- 
calculus and Year 12 calculus. Maybe I’m looking in the wrong places? I’ve got a couple of 
UK texts for the old Derive versions which I can adapt for things like applications of Calculus, 
etc. But that only gets me to week 3 in the course. 

My class is once a week, 3 hours late afternoon/early evening. I do a lecture up till 6, then a 
1 hour lab after a break. Ideally I’m looking for a couple of lab problem solving exercises 
each week that will consolidate the theory stuff after the tea break and stop them sneaking 
out of the lab when my back is turned! Have you got any ideas? I’ve looked at the Derive 
User Group discussion but it all seems so esoteric and difficult! 

I’d appreciate any feedback and suggestions as to what I could do. 

Regards 
Martin Lindsay 

If you have any advice for Martin Lindsay then please contact him. Many thanks, Josef 
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January, February and March were very busy with receiving, answering and for- 
warding DERIVE - related questions. I bothered and urged Albert Rich and Theresa 
Shelby with lots of questions. I hope that you all will find the questions and the an- 
swers interesting enough to be published in the DNL. Many problems address dif- 
ferences between DERIVE 5 and DERIVE 6. I recommend to check if you need to 
adapt functions and programs of your one which work properly in DERIVE 5 for 
their functionality in DERIVE 6. I can tell from my own experience that I sometimes 
had very unexpected results working with old DfW5-functions in DfW 6. Josef 

Hello Josef, 

I apologize for the delay in responding to the questions raised in your emails. Today I answered the 
question in your 18 January email via separate email. The following is in response to your numbered 
questions in the file questions_f or_Albert . dfw: 


Question 1 (presented by Josef) 

2 2 

*1: 14. n - 13. n < 24. n + 5-n + 9 

I want to simulate solving this inequality step by step manually; 

#1 -(14n A 2-13n) 

and I expect 0 < 10n A 2+18n+9 (as I ever found working with DERIVE 5). 

2 2 


#1: 

14 n - 13 n < 24 n + 5 n 

+ 9 

#2 : 

2 2 
(14 n - 13 n < 24 n + 5 n 

2 

+ 9) - (14 n - 13 n) 

#3 : 

2 

0 < 10 n + 18 n + 9 


#2: 

2 2 

(14 +n - 13 +n < 24 +n + 

2 

5 + n + 9) - (14 + n - 13*n) 

*3: 

9 

n + (5 + n + 9) > - 

2 


and another attempt is even worse: 


*4: 

2 2 

(14 *n - 13 *n < 24 *n + 

2 

5 *n + 9) - 14 *n + 13 *n 


2 

10*n - 164 *n + 9 

#5: n > - 

132 

Answer 1 

In Derive 6 when an expression containing inequalities is simplified, like all other functions and opera- 
tors (except equalities as discussed in Question 5 below), the inequalities are actively simplified as 
well. For example, 

2 + 3 

simplifies to the scalar value 5 and 

2 < 3 

simplifies to the Boolean truth-value "true". Therefore, 

(2 < 3) - 2 

simplifies to "true - 2", since it makes no sense to add scalars to Boolean truth-values. 
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Question 2 (presented by Jan Vermeylen) 

Why is vector (factor (k) , k, 2005 , 2006) simplifying without factoring 2006 into its prime fac- 
tors? 

*6: FACTOR (2006) = 2.17*59 

but: 

* 7 : VECTOR (FACTOR (k), k, 2005, 2006) = [5*401, 2*1003] 

or 

k [2 3] 

*3: VECTOR (FACTOR (10 ), k, 1, 3) = 1.2.5, 10 , 10 J 

#9: FACTOR ([10 , 100, 1000]) = 10*[1, 10, 100] 

Answer 2 

The problem factoring integers other than at the top-level has been resolved for the next release 
of Derive 6 as discussed in my email to the Derive Newsgroup. 

Question 3 (presented by Andre Schmidt) 

m + x 

The D£'7?/F£’-calculation J J* * dx = — 

m 

Are there any explications for this behaviour? 

Answer 3 

The antiderivative INT (#e A (x*m), x) simplifies to #e A (m*x)/m - 1/m instead of just #e A (m*x)/m be- 
cause the limit of the latter as m approaches 0 is infinite. Whereas the limit of the former as m ap- 
proaches 0 is x, which is the same as the result of substituting m = 0 before integrating the above 
expression. This is why Derive includes the constant -1/m in the antiderivative. Note that similar re- 
marks apply to the antiderivative returned by INT(x A m, x) where the problem occurs at m = -1 . 


shows a strange (incorrect?) result. 

m 


Factoring out the gcd of the 
elements of the vector is a 
nice new feature!! 


Question 4 (presented by Fritz Mohr) 


*11: M := 


0.4 0.2 

0.5 0.3 


0.174 0.094 ' 
0.235 0.127 . 

but defining n as integer and positiv 

#13: n :■€ Integer (0, «) 


*12: [N*M, 


0.26 0.14 ' 

0.35 0.19 . 


n 

*14: N 


n n 


0.4 0.2 

n n 
0.5 0.3 


This is obviously wrong! DERIVE 5 works correct! 
Do you have any advice? 


using the VECTOR-comrriand gives no problems: 


n 

' 0.26 0.14 ' 


' 0.174 

0.094 ]' 

VECTOR(W , n, 2, 3) = 

. 0.35 0.19 . 

j 

. 0.235 

0.127 J. 


Answer 4 

Even if the degree of a matrix power is symbolic, the next release of Derive will not distribute the 
power over the elements of the matrix. 
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Question 5 (again presented by Josef) 

I am co-author of a series of technology based textbooks for secondary schools and wanted to support a col- 
league in treating simultaneous linear equations. See what happens: 


#16: 

eql := 

3 * x - 2 < 

y + z = 10 

#17: 

eq2 := 

4 * x - y 

+ 2 * z = 20 

#18: 

eq3 := 

2 + x + y 

- 5*z = -10 

#19: 

[eq4 

= 4 +eql - 

- 3*eq2, eq5 := 2*eql - 3*eq3] 

#20: 

eq4 = 

(- - 

■ 2 * z = -20) 

#21: 

eq5 = 

(17. z - 

7 *y = 50) 


This is ok. But collecting the new equations in a vector - strange output!! 

7 *x - 3 *y + 3 *z = 30 
5*x-y-4*z=0 . 

This does not happen in DERIVE 5. 1 am sure that this might confuse many teachers and students. 

Answer 5 

Derive 6 artificially suppresses active simplification of equalities entered at the top-level, making it 
possible to manually solve them by adding and multiplying constants. However, this suppression 
leads can lead to confusing results like those raised in Question 5. For version 7 of Derive, some 
other method will have to be found allowing equations and inequalities to be solved manually without 
suppressing the simplification of these expressions. 


# 22 : 


eq4 

. eq5 . 


Question 6 (presented by Gary Turner, Rochester College, and others) 

I have a question concerning rectangular coordinates and spherical coordinates. 

Using my calculus text, I convert the spherical point (9, n/4, n) to the rectangular point (0,0,9). However, when I 
graph these two points using DERIVE , the points do not match up. That is, I’m seeing two points in space when I 
should be seeing only one point. I’m graphing (9, n/4, n) in the spherical system and I’m graphing (0,0,9) in the 
rectangular system. 

Any help would be greatly appreciated 

The next minor release of Derive will convert 3D data point plots using the current coordinate 
specified by the user. I apologize for the inconvenience this bug may have caused. 

Is it possible to avoid simplification of hyperbolic functions. A colleague wants to apply 

sinh(x+y) = sinh(x)cosh(y)+cosh(x)sinh(y) in both directions and expand sinh(x)-sinh(y) into a product without 

simplifying the expressions in expontial functions. 

Is this possible? Are there any tricks or settings? 

No, in the current version of Derive hyperbolic functions are always converted to exponential form 
and are not converted back to hyperbolic functions for display. However, I hope to do this con- 
version for the next major release of Derive. 

By the way, a user from Austria asked, if it is possible to programm the graphic commands (plot or draw ....) 

You know that this is one of our great wishes for the future - it is just a reminder 

Yes, the ability to program plot commands is high on the priority list for the next major release of 
Derive. 

Hope this helps, 

Albert Rich 
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Question 7 (presented by Eberhard Lehmann, Berlin) 

Dear Josef, once more the best for 2005 together with a question to the /)//A7r/:'-specialist: 
Why does it need so much time for calculating: 

dreia(n) = 

If n = 1 
19 

#23: If MOD(dreia(n - 1), 2) = 0 

dreia(n - l)/2 
3*dreia(n - I) + 1 


#24: [TABLE (dreia(n) , n, 1, 10)]' 


T 1 

2 

3 

4 

5 

6 

7 

S 

9 

10 V 

.L 19 

SS 

29 

SS 

44 

22 

11 

34 

17 

52 J. 


takes 0.05 sec 

#26: [TABLE (dreia(n) , n, 1, 22)]' 

takes 220 sec!! 

Answer 7 

Hello Josef, 

The time required to compute dreia (n) using the definition 

dreia(n) := IF(n = 1,19, IF(MOD(dreia(n - 1), 2) = 0, dreia(n - 1)/2, 3 dreia(n - 1) + 1)) 

increases dramatically with n because of what I call the "Fibonacci Fan-out Problem". 

This problem is discussed in the Derive 6 on-line help in the section on Derive programming called 
"Recursive Functions". Consider the following definition of a function for computing Fibonacci num- 
bers: 


FIB(n) := IF(n < 2, n, FIB(n - 1) + FIB(n - 2)) 

FIB is slow because to compute FIB(n), FIB(n-2) is computed twice: once for FIB(n) and once for 
FIB(n-l). Similarly, FIB(n-3) is computed three times: once for FIB(n-l) and once for each of the two 
computations of FIB(n-2). In general, FIB(n-i) must be computed the i th Fibonacci number of times. 
For example, FIB(26) re-computes FIB(1) 75025 times! Thus the number of redundant calculations, 
and hence the time required to compute the nth Fibonacci number using FIB increases rapidly (Fibo- 
naccially?) with n. 

Such redundant computation is typical for recursive definitions that naively implement recurrences, 
like Fibonacci, that are doubly recursive. In the case of the dreia function, the solution to the problem 
is to simply save the value of dreia(n-1 ) so it can be reused as follows: 

dreia(n) := PROG(IF(n = 1, RETURN 19), n := dreia(n - 1), IF(EVEN?(n), n/2, 3-n + 1)) 

Using this definition, it takes less than a 1/10 of a second to compute [TABLE(dreia(n), n, 22)]'. An 
even more efficient definition of dreia can be implemented using the LOOP programming construct. 

Hope these explanations help. 

Aloha, 

Albert D. Rich 
Co-author of Derive 
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A Hardware Problem reported by Tania Roller and resolved by Gunter Schodl: 


Tania had a problem appearing on her brandnew laptop 
computer. Working with DERIVE she faced problems 
with the screen and found screen appearances as shown 
in the right picture . 



At the occasion of a seminar she asked Gunter who is one of the most experienced IT-experts I know and he 
immediately could help. Tania's graphic card was too fast. So they set ’’Hardware acceleration" back using a 
slider bar. Start navigating through the menus with the "Display Properties" > "Advanced" > "Troubleshoot" and 
reach finally "Hardware acceleration". ( Thanks to Walter Wegscheider for providing the English "path”!) 

In German: > Eigenschaften >Einstellungen > Erweitert > Problembehandlung > Hardwarebeschleunigung 



A mail from Canada 

McDougall Duncan duncanemcdougall@hotmail.com 

Dear Josef, 

Getting back into my routine and getting things back on schedule took more time than I anticipated 
this year but things have finally settled down a bit. What I want to tell you is how impressed my asso- 
ciates are with DERIVE. They had worked with MAPLE and couldn't believe how simple and direct 
using DERIVE was. WE ARE IMPRESSED ! ! ! 

A mail from Sweden 
Dear Josef! 

Thank You for your mail. 

I've got DNL # 3 7-52 and now Tm happy to get the first 36 issues too. 

That's great! There is always something interesting in each issue so please continue to publish them 
on the internet. I'll also appreciate that You rewrits it for Derive 5. 

Best regards from 
Sture Farnstrom 
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Thierry Dana-Picard 


Dear Theresa, 

For a couple of months I have worked with Derive 6. 1 and enjoy it very much. 

I would like to share with you three cases where Derive gives the correct answer when asked to give a 
result immediately, but the step-by-step indications should perhaps be improved (IMHO): 

1. the identitites employed here are not trivial; an average student would not understand them. 

2. the same remark is valid here. In the other direction, the fact that Derive computes this integral for 
a "general" parameter r is remarkable. 

In both cases, a shorter way exists. 

3. here an explanation is needed: what leads the student to the usage of this (non trivial) formula? 

The files are attached in the corresponding order. 

May I once again thank you for the wonderful work you and your team are doing. 


Yours sincerely, 

Thierry Dana-Picard 

Dept of Applied Mathematics - Jerusalem College of Technology 


1 1 : 


f 4 * 

dx 

- + ^/(a - x) 


2 

J FC/(a + b>x) , x) dx -► SUBST x-F x, 

b U l 



dx, x, ./(a + b-x) 


*5: 


77/ 2 


0 


1 


r ■ (- LN(COS(x)) + LN(SIN(x))) 

1 + C 


dx 


b 

J (F(x) + F(a + b - x)) dx 
b a 

J F(x) dx — 

a 2 


/ 

' 2 




*’</(- x + a) 



#2: - 2. SUBST 

dx, x, 4(.-x + a) 

2 


#6: 

\ 

V(- x + a) + x ^ 




Dear Thierry, 

Thank you for your enthusiastic response to Derive 6.1. We are very glad that you enjoy 
it and appreciate your suggestions for improvement. I have forwarded these particular 
suggestions to Albert Rich, the author of the display step feature. 

It is a pleasure hearing from you. 

Regards, 

Theresa 


Hello Dr. Dana-Picard, 

I also want to thank you for your kind words and suggestions for the Derive 6 Display Step feature. 
First some general comments about this new feature: 

The transformation rules Derive displays are those it uses to simplify an expression. They may or 
may not be the same as those currently taught to students. However, if teachers see an advantage 
to an unfamiliar rule used by Derive, they may want to ask their students to verify the validity of the 
rule and then the students will have an additional tool in their arsenal. 

Thanks in part to your examples and suggestions, numerous improvements to the Display Step rules 
have been made for the next release of Derive. New algebra and calculus simplification 
rules significantly improve the understandability of the steps used to simplify your example 1 (i.e. the 
double integration of the arccosine of x). The resulting dfw file is attached to this email. 
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Your original complaint for example 1 is that "the identities employed here are not trivial; an average 
student would not understand them". Instead of re-deriving integration rules for each new problem 
(e.g. using integration by parts or substitution), Derive uses general purpose rules. Rather than being 
a mysterious "black art", integration is reduced to determining the right rule to apply based on the form 
of the integrand. 

As I suggest above, students encountering a rule with which they are not familiar provides the teacher 
with the perfect opportunity to ask the students to verify the rule. It seems to me that the ability of 
students to derive general purpose rules is preferable to their re-deriving special cases of those rules 
each time a new problem is encountered. Also the recognition that there are general purpose rules 
may be enlightening to some. 

For your examples 2 and 3 the definite integration rule in question is 


INT(F(x), x, a, b) ~> INT(F(x)+F(a+b-x), x, a, b)/2 


While apparently not yet widely known, it may be a rule that teachers should know about. As to your 
question as to when it should be applied for a given integrand F(x), simply compute F(x)+F(a+b-x) and 
see if it is simpler than F(x). Note that before simplifying F(x)+F(a+b-x), the domain of x should be 
declared to be the open interval (a,b) . In both your examples, the F(x)+F(a+b-x) simplifies to 1 which 
is certainly easier to integrate than the original integrands. 

The symbolic power r in your example 3 unfortunately results in the display of the steps required to 
convert the expression to the exponential form used internally by Derive to represent symbolic pow- 
ers. I hope to avoid the display of these unnecessary steps in a future release of Derive. 

Hope this discussion has helped. 

Aloha, 

Albert 


Dear Albert, 

Thank you very much for your detailed answer and for the dfw file; it is helpful. 

For the rule 

INT(F(x), x, a, b) -> INT(F(x)+F(a+b-x), x, a, b)/2 

I agree with you that it is important, the point was really that it is not so well known (I use such trans- 
formations in one of my papers, but a friend of mine, in another institution, told me that he would not 
hope his students to make such work). 

For the parametric example, it was a good surprise that Derive computes the integral for general r. We 
used this example as a core example in a joint paper with J. Steiner. 

I'll try to continue being in touch. 

Yours, 

Thierry. 


Dear all, 

this is a very interesting discussion. I’d like to publish it in the next DNL. I am sure that many of our 
members - especially the teachers among them - would appreciate this. It is one of the special features 
of our bulletin, that we can spread such "Insider Information" among our world wide community. 

Best regards to Hawaii and Israel, 

Josef 
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Hello Josef, 

I would be delighted if this discussion of the Derive 6 Display Step feature was republished 
(especially considering the amount of time it took to write the original email :=). However, as 
educators with experiences broader than mine, I think you and the other recipients of this 
email should give your own insights into the full potential of the display step feature. In my 
own self interest, I volunteer you, Josef, to over-see this 'discussion group'. 

To start things off, in a slightly altered state of consciousness at the beach today, I came up 
with the idea of a Display Step option for Derive 7 that would only show the antecedent (i.e. left 
side) of the transformation rule being applied. This would then be a great opportunity for a 
teacher to ask 'Ok students, what is the consequent (i.e. right side) of this rule?' When the 
class proposes a rule, the teacher asks them to prove its validity and the domain to which it 
applies. If the proposed rule turns out to be invalid, the students will not be able to prove it 
and they will have to come up with another proposed rule. 

When the correct rule is finally proved, the teacher can send a command to all (or some) cop- 
ies of Derive on the local area network to mark this rule as "available for use". Then the Dis- 
play Step feature will use and display the rule without bugging the students any more. I think 
this is a good example of Bernhard's scaffolding paradigm for teaching mathemat- 
ics. Students taught in this way would not be forced to clutter their minds memorizing rules 
(ugl). Rather, they would learn how (if they felt it necessary) to re-derive rules encountered 
when solving real-world problems in the future. 

Just a thought. 

Aloha, 

Albert 

P.S. As far as the rule 


INT(F(x), x, a, b) ~> INT(F(x)+F(a+b-x), x, a, b)/2 

not being well known: My being forced to use it to automate definite integration seems like a pretty 
good indication of its utility to me. It should be taught to students if it is not already. I know it sounds 
kind of spooky, but maybe we should listen more to what computers have to 'say' about how the Uni- 
verse works... 


Dear All, 

Thanks to all the people who related to my original e-mail. This discussion is very interesting. 

I totally agree with Albert Rich's final remark. Even in education, too many people expect from a CAS 
to make computations where they cannot or where they do not want. 

A step-by-step feature is often understood as an assistent for "a posteriori understanding"; the pro- 
posed new feature for Derive 7 could incite educators to have another usage of Derive, for an "a priori 
study". For integration, this can provide ways towards a more profound insight. 

As far as I know, the formula under consideration does not appear, neither in textbooks nor in web- 
sites. Please tell me if I'm wrong. 

Enjoy sun and beach. Here it's real winter. 

Yours, Thierry. 
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Question forwarded by a DERIVE User: 


# 1 : 


# 2 : 


#3: 


#4: 


z := 0.7 + 0.7 -l 


r n 

1000 EULER_PHI(n) • z 

APPROX Z , 30 

n=l n 

1 - z 


r n ^ 

1001 EULER_PHI(n) • z 

APPROX I , 30 

n=l n 

1 - z 


APPROX 


1002 EULER_PHI(n) • z 

I . 30 

n=l n 

1 - z 


= -1.695953134 + 0 . 005234424543 • L 


= -1.675277045 + 0 . 02591170059 - L 


= 189.6955839 - 1 . 188037185 • L 


Why does this jump so much? 
C 

#5: APPROX 


1001 EULER_PHI(n) • z 

I . 30 

11=1001 n 

1 - z 


= 0.02067608847 + 0 . 02067727604 - L 


#6: APPROX 


#7: APPROX 


1002 EULER_PHI(n)-z 

I . 30 

n=1001 n 

1 - z 

n 1 

1002 EULER_PHI(n) • Z 

I , 30 

11=1002 n 

1 - z 


= 0.02067555186 + 0. 03402483983 - L 


5.366188532-10 + 0 . 01334756378 - L 


There's nothing in these areas to justify the jump in the above #4. 
The Euler_Phi() function seems to be working fine, the problem is 
somewhere in the summation. What’s up here? 

Albert Rich answers: 


If the difference between the upper and lower limits is less than a 1000, Derive approximates definite 
sums by simple iteration. If the difference is greater than a 1000, Derive attempts to perform the 
summation by finding the antidifference of the summand (for details, see the on-line help topic "Calcu- 
lus > Sum command"). Unfortunately, in your example, this results in huge round-off errors. 

The optional fifth argument of the SUM function specifies the step size to be used for the summation. 
The default step size is 1. If an explicit step size is given, simple iteration is used for the summation 
no matter how large the difference in the limits. Therefore, you can obtain an accurate sum by includ- 
ing a step size of 1 in the call on SUM as follows: 

APPROX(SUM(EULER_PHI(n)-z A n/(1 - z A n), n, 1, 1002, 1), 30) 

which simplfies to 

-1 .6752775822038766521 1 335615807 + 0.039259264377691 1 303838052332962*#i 


We apologize for the inconvenience. 

Hope this helps. 

Aloha, 

Albert 
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Jaime Marcos 

Subject: Unexpected Complex Values in Factorial 
Dear Derivians, 

With the DfW 5.06 if you try to get, for instance, (-40.625)! in exact mode, you get a rational expression times 
(3/8)!, that simplifies in turn (working with 25 significant digits) to 1.6 669898 143291 70810 E-47. Fine. 

But if you try to get (directly) the approximate value of (-40.625)!, you get an ugly complex, 

2.8457296 1621 0276126327 153-1 0 A -47 - 1.178739801881079045325915-10 A -47 -i 

Does it happen the same thing in DfW 6.1? 

Best wishes for the 2005, 

Jaime Marcos 


Hello Jaime, 

Thank you for pointing out this defect in the implementation of the algorithm Derive uses for approxi- 
mating the factorial operator and Gamma function. Investigating the problem revealed that it only 
occurs when approximating the factorial of negative numbers greater than a certain magnitude. The 
magnitude is dependent on the digits of precision being used. For example, at 10 digits of precision, it 
occurs for numbers less than -19.6. At 25 digits of precision, it occurs for numbers less than -40.4. 

The defect has been found and corrected for the next release of Derive 6. Until then, when approxi- 
mating the factorial of negative numbers, increase the digits of precision sufficiently to avoid complex 
results. I apologize for the inconvenience. 

Aloha, 

Albert D. Rich 
Co-author of Derive 


We have a new member from Spain who started programming with DERIVE. He worked through the 
JOSEPHUS Problem (first presented by Riideger Baumann) from earlier Newsletters (DNL#52, #53) 
and found a QUICKER-JOSEPHUS-Solution. You can compare calculation times. Angel sent also a 
small program dealing with the famous QUEENS Problem on a checker board (next DNL). 

Angel Aguirre 

Josephus problem: my solution (borrowing some ideas from Welke) 


Qui ckerJosephus(n , s, p_, 1_, u_, v_) := 
Prog 

[p_ := s, 1_ := n, u_ := [1, . . . , n] ] 
Loop 

V_ := [p_, p_ + s 1_] 

u_ := APPEND(DELETE(u_, v_) , u_|v_) 
p_ := s - M0D(1_ - p_, s) 

1_ DIM(v_) 

If 1_ = 0 
RETURN u_ 
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Finally we have a Tl-related question: 

Lieber Josef, 

Einige Fragen zum TI: 

• Kann man vom Homefenster auf eine Zelle im Data-Matrix-Fenster irgendeiner Datei zugrei- 
fen? 

• Wenn ich in einer Spalte einer Data-Matrix Funktionen stehen habe, kann ich die dann elegant 
plotten? 


Dear Josef, 

• is it possible to address one single cell of a DATA-Matrix in the HOME-Screen? 

• having functions in a column of a DATA-Matrix, is there an easy way to plot the func- 
tion graphs? 


I create a data-sheet heiner (what else) and fill the The columns can be transferred into the home 
first two columns with functions and their deriva- screen as lists. 


tives 


PATA 

1 I I 1 





cl 

c2 

c3 

g4 

g5 


1 

X 

1 





2 

2*x^2 

4+x 





3 

3*x^3 

9*x^2 





4 


16*x A 3 





5 

5*x^5 

25+x'M 





6 

G+x^G 

36*x A 5 





7 

7*x^7 

49+x^G 





c2= 

d< cl. 

X> 


MAIN 


RAD EXACT 

FUNC 1 


Single elements can be selected using a second 
index (in brackets!). 


itiTTiT fs^ y f::t y fht y fe Y fa^ y "i 
I t fr— |fl 1 gebra |Ca 1 g lother |Prgn 1 0 |c 1 ean Up | 

■ heinert2] 

{1 4 ■ x 9x 2 16 ■ x 3 25 ■ x 4 

36 x 5 ► 

■ (he iner[ 2]) [3] 

9 ■ x 2 

■ Graph (heiner[2])[3] 

Error: Undefined 

variable 

■ (heiner[2])[3] + yl(x) 

Done 

heiner[2] [3]+yl<x> 


MAIN RAP EXACT FUNC H/30 



Unfortunately it is not possible to plot the graphs 
directly. Defining a function in the Y=-Editor 
doesn't help. 

The function is recognized, because yl (2) re- 
turns the correct value 36 in the home screen. 


This is the error message after switching to the 
GRAPH-Window. 


fFiTHJY F£t Y F3 Y FH Y FE^ Y FfiT YfF 

| t f— IZoon |T raoe IRegraph |Math |Draw | t 

sm 


i 

r 


f ERROR ■> 



Undefined variable 







^CEnt.er=G0T0!> 

(!ESC=CRHCEO ; 





main 

RAP EXACT 

FUNC 



You must explicitly call the element and then as- 
sign it to a function. Now it works! 


I'F i"7 T TlT fet y f::t y fh^ Y FE Y FhT T -I 

| ■* |r 1 gebra |Ca 1 g |ot-her |Prgn 1 0 |c 1 ean Up | 

■ he inert 2] 

Cl 4 ■ x 9 ■ x 2 

■ (he inert 2]) [3] 

■ 9x 2 + yl(x) 

16 ■ x 3 25 ■ x 4 36x 5 ► 

9 ■ x 2 
Done 

|ans<l>+yl<x> 

i i 


I MAIN RAP EXACT FUHC 5/50 


fFiTHQY F2w T FS T FH Y FEt Y FfiT YF? jf Y r "1 

K Zoom |Traoe Regraph |Mat-h|Draw|T ^ |:( | 







MAIN RAP EXACT FUNC 
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GEOMETRIE MIT TI-92 UND TURTLE-GRAPHIK 
Geometry with TI-92 and Turtle Graphic 

Mandala von Westenholz, Osnabrtick 

Im Anschluss an die Ergebnisse der TIMSS (Third International Mathematics and Science Study) 
und der PISA-Studie hat das Zentrum fur Kreativitatsforderung in Osnabriick mit jungen Schii- 
lern eine Untersuchung gemacht, NEUES LERNEN in Mathematik zu erproben. Hier sind einige 
Programme, die ich als jiingste Teilnehmerin an diesem Projekt (ich bin 1 1 Jahre alt und besuche 
ab Herbst 2003 die gymnasiale Klasse 8), erstellt habe und die illustrieren, wie mit Hilfe der 
LOGO-Schildkrotengeometrie (TI-92PLUS - Turtle Grafik [1] ) das Grundverstandnis fur die Geo- 
metric sowie die geometrische Anschauung entwickelt und gefordert werden konnen. 

After presentation of the results of TIMSS and PISA the Creativity Support Center in Osnabruck 
started an investigation to test NEW LEARNING in mathematics. Here are a few programs which I - 
as the youngest participant in this project (I am 11 years old and will attend the 8 th form in a gymna- 
sium next fall) produced and which shall demonstrate how basic understanding of geometry can be 
developed and supported by the use of LOGO-Turtle geometry). 

Hauptdarstellerin ist die hier abgebildete Schildkrote TURTLE- 
FIX. Unter ihrer Mitwirkung sollen die folgenden Probleme 
gelost werden: 

Turtle TURTLEFIX is our main actor and it will help solving the 
following problems. 

Problem 1: Definiere eine Prozedur „square“, die ein Quadrat der Seitenlange x erzeugt. Wie ver- 

wendet man dieses Programm, um eine Familie von „Wachstumsquadraten“ zu erzeu- 
gen? 

Problem 1: Define a procedure “square", which generates a square with side x. How to use this 

program to create a family of “Growth-Squares"? 



■ FiT Y Y FjT Y FHt YfFSy FfiT 

Misc|Lessons PentT|HidezShow|tz |Polygon 

i ► 

■ tg\ ini ttrt 1 ( " 1 : 1" ) 

Done 

■ Q + ^turtle G.0G 

■ O + xturtle 0.0G 

■ O + yturtle O. 00 

■ square(40) Done 

i 

i 


tsr\ShowTrtlO 

TQ DEG HF'F'RDK PfiR E/HD 



rriT^nT ^ fe y f6tY 

^ir— IContnol lzo|Uan|Find...|Mode| 

: quadrate^x} 

: Pngn 

: If x > 1 0 Then 
: squane<x> 

: quadrate <x-5> 
: Endlf 
: EndPrgn 


TG DEG RPPRDK 

PAR 



Die abgebildete Familie von Wachstumsquadraten habe ich mit quadrate (80) aufgerufen. 
I produced the family of Growth-Squares by executing quadrate (80) . 
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(Before leaving the turtle alone on the screen One has to set the initial values for the little creature: 
xturtl e and yturtl e are the coordinates of the starting point and 0turtl e is the initial direction, 
Josef) 

Diskussion: Das entsprechende Programm in der Programmiersprache LOGO lautet: 

This is the respective program written in LOGO: 

PR Quadrate :x 
wenn :x < 10 dann rk 
Quadrat : x 
Quadrate :x - 5 
ENDE 

PR Quadrat :x 
wh [ vw :x re 90 ] 

ENDE 

(mit den Schildkrotenbefehlen wh = wiederhole, vw = vorwarts, re = rechts, x = Seitenlange) 

(Turtle commands wh = repeat, vw = forwards, re = turn right, x = side length) 

Als nachstes stellt sich die Frage: Wie kann ich mit dem TI-92 und Turtlefix ein regelmaBiges Fiinfeck 
konstruieren? Mit der Programmiersprache LOGO untersuche ich hierzu, wie Turtlefix fur eine ent- 
sprechende Schildkrotenrundreise zu programmieren ist. 

Next question is how to guide Turtlefix to plot a regular pentagon on the TI-92? At first I used LOGO to 
investigate how to program Turtlefix for a respective turtle roundtrip. 

Die Figur zeigt: Die Schildkrote dreht iiber die Aufienwinkel. Die Winkel, die wir normalerweise mit 
dem Winkelmesser (Geodreieck) messen, sind Innenwinkel. Insgesamt muss sich Turtlefix bei einem 
Umlauf um 360 Grad drehen, also bei jeder Ecke um 

p = = 72 ° fur ft = 5. 

ft 

Ein AuBenwinkel ergibt sich somit als Quotient aus 360° durch die Eckenanzahl n. 

The figure shows: the turtle turns according to the 
outer angle (left picture). Angles measured by our 
tool are usually inner angles. Turtlefix must turn 
about 360 degrees making one complete penta- 
gon, so at each vertex turn about 36075 = 72°. 

The right picture shows the wrong programmed 
Turtlefix. 



Ich erzeuge die Turtleflx-Pentagonflgur mit einem iterativen LOGO-Programm (iterativ = sich wie- 
derholend) dann so: 

PR PENTAGON :x :n 

wh :n [ vw :x re 3 60/n ] This is the LOGO program. 

ENDE 

Problem 2: Wie ist mit dem TI-92 und der Turtle-Graflk ein interaktives Programm zu erzeugen, 

dass Turtlefix ein Pentagon zeichnet? 

Problem 2: How to design a program with the TI-92 and Turtle-Graphics to make Turtlefix drawing a 

pentagon? 
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Losung / Solution: 


T I 

F£t 

Control 

ixo 

fFH^Y 

Uar 1 

F£ T F6t Y 

- ind... Mode 

: pentagontx) 
: Prgm 

: For i , 1 , 5, 1 
: forwardtx) 
:rt<72> 

: EndFor 
:DelUar i 
: EndPrgml 


TG DEG APPRDH PAR 


iTi^mr FE ' ¥ ' T F5 T FH T F£t 

t f— - Zoom Trace Regraph Math 

FfiT X- Y =:! 

DrauH k 

“ ■■ ■-** •-■■-■II JL ■=• 

■ pentagon(50) Done 

■ showtrtlO Done 

■ "30 + xturtle : 0 > 

G 

■ 9G * Sturt le 9Q 

■ pentagon(30) Done 

■ hidetrtlO Done 


rr\ 

0 

^ / 



TG DEG APPRDK PAR 


Beispiel: Das Netz eines Pentagon-Dodekaeders lasst sich nun programmieren (Turtlefix ist ver- 
steckt!). 

Example: The net of a pentagon dodekaeder can be programmed. (Turtlefix is hidden!) 



This net is not produced by a program. As you can 
see it is made using Cabri on the TI-92. 

Mandala did it this way too. It is up to you to write a 
program using pentagon ( ) . 


Aus dem regelmaBigen Fiinfeck ensteht durch Verlangerung seiner Seiten das Pentagramm. Diese 
Figur war schon im Griechenland der Antike von Bedeutung, insbesondere als Symbol der Bruder- 
schaft der Pythogoreer, also der Anhanger des beruhmten Mathematikers Pythagoras. 

Lengthening the sides we obtain a pentagram. This figure had its importance in antique Greece as 
symbol of brotherhood of Pythagoreans. 

links: Griechische Miinze mit eingepragtem Pentagramm (420 
v.Chr) 

rechts: Griechische Miinze mit Portrat von Pythagoras (430 v. Chr.) 


left: Greek coin with pentagram (420 BC) 
right: Greek coin with picture of Pythagoras (430 BC) 

Problem 3: Ein Archaologe findet eine antike griechische Miinze. Das eingepragte Pentagramm ist 

beschadigt. Um die Figur dennoch im Internet archivieren zu konnen, programmiert er 
das Pentagramm mit dem TI-92 und Turtlefix. Wie muss das Programm beschaffen 
sein? 

Problem 3: An archaeologist finds an antique greek coin showing a damaged pentagram. He wants 

to archive the figure in Internet and programs the pentagram on his TI-92 supported by 
Turtlefix. How will his program look like? 




Losung / Solution: 


The instruction how to draw a pentagram leads to a program (fixed side length 60), star(-30,-20,60) is 


slightly improved. 



star () 

Prgm 

For i ,1 ,5,1 
forward (60) 
rt (144) 
EndFor 
EndPrgm 


fFiT^T F £ T T F5 T FH Y FEt T FfiT TT? x>Y=:i 

t f— |Zoom |T race IRegraph |Math |Draw | ■* ^ |i! 

■ showtrtlO Done 

■ start 0, -40, 40) 

Done 

■ start! "3D, -20,60) 

Done 


< / "*■ 
l/ 



TG DEG APPRDX PAR 


D-N-L#57 


Mandala von Westenholz: Turtle Graphic on the TI-92 


p45 


Fur einen Ausflug in die Welt der fraktalen Geometrie von Benoit Mandelbrot mache ich von re- 
kursiven Programmen Gebrauch. Was versteht man hierunter? 

„Rekursiv“ = „zurucklaufend“, das bedeutet: Eine geometrische Konstruktion oder eine bestimmte 
Beweisfuhrung wird verschiedene Male wiederholt. 

Problem 4: Wie ist eine rekursive Prozedur beschaffen (eine Prozedur, die sich selbst aufruft), die 

das Fraktal eines binaren Baums erzeugt? 

Problem 4: Which is the recursive procedure (a procedure which calls itself), which creates the 

fractal of a binary tree? 


Losung / Solution: 

Fur den Binarbaum gilt: ein Teil ist eine exakte Kopie. 
Each part of the binary tree is an exact copy of the whole. 


Programmiere den Stamm (der Lange x) 
Programmiere den linken Teil mit halber Lange 
Programmiere den rechten Teil mit halber Lange 
Programmiere die Ruckkehr zum Ausgangspunkt 


Program the trunk (length x) 

Program the left part (half length) 
Program the right part (half length) 
Programm the return to the starting point 


fFiT^OT Y FS Y FH Y T FfiT V 7 

^ t—IZoon Trace IRegraph Math IDrauK r M 

■ G + xturtle : "50 > 

-5G 

■ 0 + Sturt le G 

■ tree(50, 45, 5) 

□one 


* 


TG DEG fiFFRDK FfiR 



Die Figur habe ich mit t r ee ( 50 , 45 , 5 ) aufgerufen. Mit t bezeichne ich die Verzweigungstiefe und 
d ist der Winkel unter dem ein Teilbaum verzweigt. 

PR TREE : x :t 

wenn :t = 0 dann vw :x rw :x rk 
vw :x re 45 
TREE : x/2 :t-l 
li 90 

TREE : x/2 :t-l 
re 45 rw :x 
ENDE 



The trees of levels 0-4 were created by using the DERIVE file logodfw5.mth which has to be adapted 
for use with DERIVE 6 and is now logodfw6.mth. It was produced as a coproduction of Josef Lechner, 
Eugenio Roanes and Johann Wiesenbauer (DNL#25 and DNL#38). Josef 
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Problem 5: Wie ist ein Koch-Fraktal mittels der Turtle-Geometrie zu programmieren? 

Problem 5: How to program a Koch-Fractal applying Turtle-Geometry? 


Losung / Solution: 

Gegeben ist eine Strecke der Lange x. Ersetze diese durch folgende Figur: 
Given is a segment with length x which is substituted by the following figure: 


• Strecke x/3, drehe um 60° nach links, 

• Strecke x/3, drehe um 120° nach rechts, 

• Strecke x/3, drehe um 60° nach links, 

• Strecke x/3. 


• Segment x/3, left turn by 60°, 

• Segment x/3, right turn by 120°. 

• Segment x/3, left turn by 60°, 

• Segment x/3. 


PR KOCH : x :t 

wenn :t = 0 dann vw :x rk 

KOCH : t-1 : x/3 

li 60 

KOCH : t-1 : x/3 
re 120 

KOCH : t-1 : x/3 
li 60 

KOCH : t-1 : x/3 
ENDE 




koch£x, y!> 

Prgn 

If y>0 Then 

kochCx/3,y-l>: lftCSO) 
kochCx/3 ? y-l>:rta2G> 
koch<xx3,y-i:>: lftCSO) 
koch<x/3, y-1 > 

Else 

forward(x> 

End If 
EndPrgpJ 


deg hfff;d:: 




"4Q ■+ xturtle 


- 10 + yturtle 


' tstartC -40, -10, 90) 

' koc-KSQ, Q) 


koch<80,0> 


go 4 etu> 
9G. 
Done 
Done 


To make it easier setting the initial conditions of Turtlefix I defined a tiny function t start (x , y , 0). 






1 -40 ■+ xturtle : -10-* yturtle : 9G ■+ 9tu^ 

9G. 

' tstartC - 4G , - 1 G , 90) Done 

■kocKSOjO) Done 

' tstartC - 4G , - 1 G , 90) Done 

1 koch(S0, 1) Done 


koch<80,l> 





DEG flPPRDK 



41 

■ tstartc -40, -10, 90) 


9U. 

Done 

■ kochCSG, G) 


Done 

■ tstartC -4G, -10, 90) 


Done 

■ kochCSG, 1) 


Done 

■ tstartC -4G, -10, 90) 


Done 

■ kochCSG, 3) 


Done 

■ tstartC -4G, -10,90) 


Done 

koch<8G.5> 

TG DEG APPRDK 

FUNC 



This is the fractal Snowflake: 



- ^^Iconiro 1 f I^oluarfF ind.. Jllodel 

:flake(x,y) 
: Prgn 

: kocIVx, y) 

: rtO20) 

: kocIVx, y) 

: rtO20) 

: kocIVx, y) 

: rtO20) 

: EndPrgn 


TG DEG APPRDK 

FUNC 



Problem 6: Ein Sierpinski-Dreieck ist ein Fraktal, das sich durch Ahnlichkeitsabbildungen (Ge- 

ometric der gymnasialen Klasse 9) so beschreiben lasst: Gegeben sei ein Dreieck ABC 
mit den Ecken (0,0), (2,0) und (1,1). Ein Sierpinski-Sieb entsteht dann durch die Defini- 
tion von drei Ahnlichkeitsabbildungen, zentrische Streckungen mit dem Streckfaktor 
k = 1/2, die, verkniipft mit Translationen um die Strecke 1 und 1/2 in x- bzw. y- 
Richtung wie folgt gegeben sind. 




D-N-L#57 


Mandala von Westenholz: Turtle Graphic on the TI-92 


p47 


Problem 6: A Sierpinski-Triangle is a fractal which can be described by similarity mappings (ge- 

ometry of form 9 in gymnasium) as follows. Given is a triangle ABC with vertices (1,0), 
(2,0) and (1,1). A Sierpinski-Sieve is generated by defining three similarity mappings, 
centric stretching with stretch factor = 1/2 combined with translations by 1 in x- and in- 
direction. 

Jeder Punkt (x,y) des Dreiecks wird in einen Punkt (x/2,y/2) abgebildet. Dadurch entsteht das Dreieck 
#1 von der halben GroBe in der linken unteren Ecke. In Matrixform geschrieben: 

Each Point (x,y) of the triangle is mapped into a point (x/2,y/2) giving triangle #1 left on the bottom 
which has half size. This can be written in matrix form. 


/: 






-» 




<y'j 




f°) 


+ 


UJ 


v^y 



(efi) Cue) (V) 


Jeder Punkt (.v.r) des Dreiecks wird in einen Punkt (x/2+l,y/2) transformiert. Es entsteht das auf die 
Halfte verkleinerte Abbild #2 des Dreiecks. In Matrixform geschrieben: 

Each Point ( x,y ) of the triangle is mapped into a point (x/2+1y/2) giving triangle #2 which has half size 
again. This can be written in matrix form: 


g- 




( x ^ 


-> 


VTy 


<y'j 


fx) 


rn 


+ 


Uy 


7y 


Jeder Punkt (x,y) des Dreiecks wird in einen Punkt (x/2+l/2,y/2+l/2) transformiert. Durch Halbieren 
aller Koordinaten und anschlieBender Translation um 1/2 entsteht die auf die halbe GroBe verkleinerte 
Kopie #3 des Dreiecks. In Matrixform geschrieben: 


Each Point (x,y) of the triangle is mapped into a point (x/2+1/2,y/2+1/2) giving triangle #3 which is a 
half size copy of the initial triangle. This can be written in matrix form: 


g- 




\y) 




\y ) 


0 


1 

27 


\y) 


+ 


(£\ 

2 

I 

V2 7 


Der iterative Konstruktionsalgorithmus lautet: Verbinde die Mittelpunkte der Dreiecksseiten, ent- 
feme das mittlere Dreieck #4, verfahre in gleicher Weise mit den iibrigen Teildreiecken #1, #2 und #3. 

The iterative algorithm is: Connect the midpoints of all sides of the triangle, remove triangle #4, pro- 
ceed in the same way with the remaining triangles #1 , #2 and #3. 

How to program a Sierpinski Triangle applying Turtle Geometry on the TI-92? 

Losung / Solution: 


Mit dem Aufruf si erpi n (50) erhalte ich den folgenden Screenshot: 
Calling si erpi n (50) I obtain the following screenshot: 


fFiTmjT y f: t fh 

t f— |Zoon |T race |Regr a 


4 

xnin= ~ 47.6 
xnax=47. 6 
xscl=l . 
ynin= - 4G. S 
ynax=4G. S 
yscl=l. 
xres=2. 

TG DEG APPRDK FUNC 




: sierpin(x) 

: Prgn 
s If x>lG Then 

: forward (x>:lft(12G}:si er p i n ( xx2 > 
: f or war d (x>:lft(12G}:si er p i n ( xx2 > 
: f or war d (x>:lft(12G}:si er p i n ( xx2 > 
: End If 
: EndPrgml 
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Das letzte Problem befasst sich mit dem 

BAUM DER ERKENNTNIS VON PYTHAGORAS 

in der Welt der Turtle-Geometrie. Jedoch: Wie sieht ein solcher Pythagorasbaum iiberhaupt aus? 
Zunachst betrachte ich vier verschiedene „Momentaufnahmen”, die Turtlefix beim Zeichnen einer 
Pythagorasfigur zeigen. Man erkennt die Quadrate liber den Katheten und der Hypothenuse. 

The last problem deals with the TREE OF KNOWLEDGE OF PYTHAGORAS in turtle geometry 
world. But how does this tree look like? At first I inspect four different snapshots which show Turtlefix 
drawing the figure. One recognizes squares above the sides of a right triangle. 



QUADRAT : x QUADRAT :x*cos(30) 
vw : x 
li 30 



re 90 QUADRAT :x*cos(30) 

vw : x*cos ( 30 ) 


Das einfache LOGO-Programm zeigt, wie mit Turtlefix der Pythagorasbaum zu programmieren ist: 

LERNE PYTHAGORAS : x 
QUADRAT : x 

WENN : x<l DANN re 90 vw :x rk 
vw : x 1 i 30 
PYTHAGORAS :x*cos (30) 

PYTHAGORAS :x*sin (30) 
re 30 vw :x li 90 
ENDE 

Problem 7: Welches TI-92-Programm erzeugt einen Pythagorasbaum? 

Problem 7: Which TI-92 program creates a Tree of Pythagoras? 

Losung/ Solution: 


PI jcontro 1 1 1 ?oluarlF ind.. Jriodel 1 

Pi=feE5,l 1 

UpI 1 


pytha(x) 

Frgm 

If x<2 Then 

rt<9Q> : forwarcKx} : Return 
End If 
squared x> 

forwardCx): lftC3Q> 

py tha < x*cos < 3Q > > 

py tha < x*s i n £ 30 > > 

ntC3Q> : f onward < x : lftC9Q> 

EndPrgn 

xnin= "19Q. 4 
xnax=190. 4 
xscl=l . 
ynin= "SI . 6 
ynax=81 . 6 
yscl=l . 
xne s=2.| 

■ tstarU20, -70, 0) Done 

■ pytha(40) Done 

pvtha<40> 

TG DEG APPRUK FUNC 

TG DEG APPRUK FUNC 

TG DEG APPRUK FUNC 
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Mandalas paper is based on the program package Turtle-Graphic from [1] and [6]. I try to get permis- 
sion from the authors for including this package to the files which can be downloaded woth DNL#51 , 
Josef. 


Online-Tutorials fiir Voyage 200 und DERIVE 6 

Haben Sie eine Einstellung vergessen? Fallt Ihnen Befehl nicht gleich ein, dann besuchen Sie die Tu- 
torials fur Voyage 200 (auch fur TI-92, TI-92PLUS und TI-89 einsetzbar) und das neue Tutorial fur 
DERIVE 6. Beide wurden von Walter Wegscheider zusammengestellt und werden weiter ausgebaut. 
Vorerst gibt es diese wertvollen Hilfen nur in Deutsch. 

http : //www . austromath . at/daten/ derive/ 
http : //www . austromath . at/daten/voyage200/ 


Beide Tutorials sind auch hervorragend fur den Gebrauch durch Schuler geeignet. 


Until now time both Tutorials are available only in German. 
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Take care when providing selfmade UTILITY-files 


Dear Josef 

I have been trying to load mth. Files as utility-files. However in some cases it works perfectly, in other cases, I 
cannot use the functions defined in the utility file. Are there any requirements to the content of a utility-file? I 
attach two files, one which can be loaded, the other can’t. 

Best wishes 
Ove 


DNU 

Many thanks Ove for this interesting question. I am sure that this might be of interest for many mem- 
bers who like to prepare utility-files. 

I inspected the mysterious file and indeed I found the message that file Funktionsun- 
dersogelseff ,t,n) was loaded as utility file, but there was no single function present!! There were no 
problems to load the MTH-file as a MATH-file. 

/I : LOAD(C: \DOKU5\DNL\DNL57\funktionsundersffigelse.mth) 

As you might know MTH-files are saved in pure ASCII-Code but DERIVE 6 works with Uni- 
code, which is an very much extended code to also include the many mathematical charac- 
ters and the characters belonging to various languages. Look at the following function defini- 
tion: 

Skra(f, t, n) := 

Prog 

t := NUHIERATOR(f) 
n := DENOPIINATOR(f) 

If POLY_DEGREE(t, x) = POLY_PEGREE(n, x) + 1 
["Skra asymptote y="; QUOTIENT(t, n)] 

"Ingen skra asymptote" 

And then look at the same function in any textprocessor or any editor: 

SkrD00e5 (f, t,n) : =PROG (t : =NUMERATOR (f ) , n : =DENOMINATOR (f ) , IF (POLY_DEGR 
EE (t, x) =POLY_DEGREE (n, x) +1, [ [ "SkrD00e5 asymptote y="], 

[QUOTIENT (t, n) ]] , "Ingen skrD00e5 asymptote")) 

As you can see, the Danish “a” looks very strange in ASCII-Code. So I renamed the file as 
Funktionsundersoegelse .mth and Skra(f,t,n) as Skra(f ,t,n) and did the same with all 
the other “Danish” function names. The same problems occur in German using the letters A, 
a, O, 6, U, ti or IJ. The problems appear giving a "regular" filename and "irregular" function 
names. 


Dear Josef 
It works! 

I will have to be more careful using all the possibilities of the Danish language (as you probably will 
using German!) 

Still not using files as utility files there seems to be no problems! 

Thanks for your quick answer. 

Best regards 
Ove 


